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ABSTRACT 

We present combined interferometer and single dish telescope data of NH3 (J, K) = 
(1,1) and (2,2) emission towards the clustered star forming Ophiuchus B, C and F Cores 
at high spatial resolution (~ 1200 AU) using the Australia Telescope Compact Array, 
the Very Large Array, and the Green Bank Telescope. While the large scale features 
of the NH3 (1,1) integrated intensity appear similar to 850 /um continuum emission 
maps of the Cores, on 15" (1800 AU) scales we find significant discrepancies between 
the dense gas tracers in Oph B, but good correspondence in Oph C and F. Using the 
CLUMPFIND structure identifying algorithm, we identify 15 NH3 clumps in Oph B, and 
3 each in Oph C and F. Only five of the Oph B NH3 clumps are coincident within 30" 
(3600 AU) of a submillimeter clump. We find vlsr varies little across any of the Cores, 
and additionally varies by only ~ 1.5kms~^ between them. The observed NH3 line 
widths within the Oph B and F Cores are generally large and often mildly supersonic, 
while Oph C is characterized by narrow line widths which decrease to nearly thermal 
values. We find several regions of localized narrow line emission {Av < 0.4kms~^), 
some of which are associated with NH3 clumps. We derive the kinetic temperatures 
of the gas, and find they are remarkably constant across Oph B and F, with a warmer 
mean value (T^ = 15K) than typically found in isolated regions and consistent with 
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previous results in clustered regions. Oph C, however, has a mean Tk = 12 K, decreas- 
ing to a minimum Tk = 9.4 K towards the submillimeter continuum peak, similar to 
previous studies of isolated starless cores. There is no significant difference in temper- 
ature towards protostars embedded in the Cores. NH3 column densities, A^(NH3), and 
abundances, ^(NHs), are similar to previous work in other nearby molecular clouds. 
We find evidence for a decrease in ^(NHs) with increasing N(H.2) in Oph B2 and C, 
suggesting the NH3 emission may not be tracing well the densest core gas. 

Subject headings: ISM: molecules - stars: formation - ISM: kinematics and dynamics - 
ISM: structure - radio lines: ISM 



1. Introduction 



Stars form out of the gravitational collapse of centrally condensed clump^ of dense molecular 
gas. Recent years have seen leaps forward in our understanding of the structure and evolution 
of isolated, star form ing clumps. Most star formation, however, occurs in clustered environments 
(jLada &: Ladall2003l ). These regions are more complex, with complicated observed geometries, and 
contain c lumps which tend to have higher densities and more compact sizes than those found in 



isolation (Motte et al. 



199S 



Ward-Thompson et al.l 120071 ). It is likely that due to these differences 



the star formation process in clustered regions proceeds differently than in the isolated cases. 
Characterizing the physical and chemical structures of these more complicated regions are thus the 
first steps towards a better understanding of the process of clustered star formation. 

It is now clear that molecular clumps become extremely chemically differentiated, as many 
molecules commonly used for tracing molecular gas, such as CO, become severely de p leted in the 
innermost regions through adsorption onto dust grains [see, e.g., iDi Francesco et al.l (|2007l ) for a 
review]. An excellent probe, therefore, of dense clump interiors is the ammonia molecule (NH3), 
with a relatively high critical density {ricr ~ lO^cm"^ for the (1,1) and (2,2) inversion transitions) 
and apparent resistance to depletion until extreme densities arid low temperatures are reached in a 
starless core's evolution (jTafalla et al.ll2004l : lAikawa et al.ll2005l : iFlower et al.ll200a ). The additional 
kinematic information provided by line observations are complementary to continuum observations 
of emission from cold dust, and the ammonia molecule in particular allows the determination of the 
gas kinetic ternperat ure and density structure due to hyperfine transitions of its metastable states 
(|Ho k Town^ll983l ). 



The nearby Ophiuchus molecular cloud, containing the dark L1688 region, is our closest 
example of ongoing, clustered star formation. The central Ophiuchus cloud has been surveyed 



^In this paper, we call prestellar objects 'clumps' instead of 'cores' to avoid confusion with the Ophiuchus 'Cores' 
discussed here. 
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extensiv ely in millimeter (I Young et alj l2006l : IStanke et al.l 120061 : iMotte et al.l Il998l ) and submil- 
limeter ([Johnstone et al.1 12004| . l2000bl ) continuum emission. These observations have revealed a 
highly fragmented complex of star forming clumps with masses M ~ 0.2 — 6Mq, the majority 
of which are embedded wi thin larger, highly fragmented structures, called 'Cores' for historical 



reasons ( Loren et al 



tion (Johnstone et al 



199d) and named A through I, which resi de only in areas of high extinc- 



2004 



Young et al. 



2006 



(sub)millimeter clum ps, ^ 40 — 50 Mr? 



the molecular cloud (I Young et al 



2006 



Enoch et al.l 120071 ). The total mass of the distinct 
makes up only 0.5 — 3 % of the total ~ 2000 M0 mass of 



Johnstone et al.ll2004l ). 



Most recent estimates put the distance to the central L1688 cloud region (often also called '/) 



Oph') at 120 pc (ILqinard et al 



2008 



Lombardi et al.ll2008l : iKnude &: Hog]|l998l ). in agree ment with 



some older results (jde Geus et al.lll989l ). but a clear consensus has not yet been reached. iMamajek 
(j2008l ). for example, find a di stance of 135 p c towa rds the cloud using Hipparchos parallax data, 
while VLBA observations by iLoinard et al.l (|2008l ) suggest that the Ophiuchus B core may be 
further distant than the rest of the cloud, at 165 pc (we also note that Oph B consists of three sub- 
Cores, Bl, B2 and B3, d escribed further i n §3). This distance is outside the range in median cloud 
thi ckness deter mined bv lLombardi et al.l ( 20081 ) of 28^'^^g pc, but in agreement with an older result 
bv lChinil (119811). The stars u sed to determine the distance to Oph B may, however, be background 
stars (jLombardi et al.ll2008l ). In the following, we adopt the 120 pc distance to the entire central 
Ophiuchus region. 

In this work, we discuss the results of high resolution observations of NII3 (1,1) and (2,2) in 
the Ophiuchus B, C and F Cores to study the distribution, kinematics and abundance patterns 
of the Cores and associated embedded clumps. We find that although the Cores are embedded in 
the same physical environment, they present very different physical characteristics. We discuss the 
observations and the combination of interferometer and single dish data in §2. In §3, we present 
the data, and detail the hyperfine line fitting procedure and derivations of kinetic temperature Tk, 
NH3 column A^(NH3) and space density n(H2) in §4 (see also Appendix A). We discuss the results 
in §5 and summarize our findings in §6. 



2. Observations and Data Reduction 



Figure [T] shows the central Ophiuchus region in 850 fiui continuum emission first mapped with 
the Submill imetre Common User Bo lometer Array (SCUBA) at the James Clerk Maxwell Telescope 
(JCMT) by I Johnstone et al.l (l2000bl) and rec e ntly r e-reduced and combined wit h all other SCUBA 
archive data in the region bv lJorgensen et al.l (|2008l ). following the description in lKirk et al.l (|2006l ). 
The Oph B, C and F Cores are labelled, and boxes show the approximate areas we mapped at the 
Green Bank Telescope (GBT), the Australia Telescope Compact Array (ATCA) and the Very Large 
Array (VLA). The details of all astronomical observations are described below. Table [1] lists the 
lines observed and their rest frequencies. 
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2.1. Green Bank Telescope 

Single-dish observations of emission from the NH3 (J, K) = (1,1) and (2,2) inversion lines, C2S 
Jat = 2i — Iq and HC5N J = 9 — 8 in the Ophiuchus Cores were obtained using the 100 m Robert C. 
Byrd Green Bank Telescope (GBT), located near Green Bank, WV, USA. The observations were 
done in frequency-switching mode, using the GBT K-band (upper) receiver as the front end, and 
the GBT spectrometer as the back end. This setup allowed the simultaneous observation of all 
lines in four 50MHz-wide IFs, each with 8192 spectral channels, giving a frequency resolution of 
6.104 kHz, or 0.077 km s"^ at 23.694 GHz. 

The data were taken using the GBT's On-The-Fly (OTF) mapping mode, using in-band fre- 
quency switching with a throw of 4 MHz. In OTF mode, a map is created by having the telescope 
scan across the target in Right Ascension (R.A.) at a fixed Declination (Dec), or in Dec. at a 
fixed R.A., writing data at a predetermined integration interval. The maps of Oph Bl and B2 
were made while scanning only in R.A. at a fixed Dec, while for subsequent targets (Oph B3, C 
and F) the scanning mode was alternated to avoid artificial striping in the final data cubes. No 
striping, however, is apparent in the final Bl or B2 images. At the observing frequency of 23 GHz, 
the telescope beam was approximately 32" FWHM. Subsequent rows or columns were spaced by 
13" in Dec. or R.A. to ensure Nyquist sampling. Scan times were determined to ensure either 
one or two full maps of the observed region could be made between pointing observations. For all 
observations, pointing updates were performed on the point source calibrator 1622-254 every 45 - 
60 minutes, with corrections approximately 2 — 3". The average telescope aperture efficiency rj^ 
and main beam efficiency rj^b were 0.59 it 0.05 and 0.78 it 0.06 respectively, determined through 
observations of 3C286 at the start of each shift. The absolute flux accuracy is thus ~ 8%. The 
average elevation of Ophiuchus for all observations was approximately 26 degrees. 

System temperatures (Tsys) varied between 48 K and 92 K over the observation dates with an 
average Tgyg ~ 62K. Table [2] gives the area mapped in each region and the final rms sensitivity in 
K (Tmb). 

Initial data reduction and calibration were done using the GBTIDlE package. Zenith opacity 
values for each night were obtained using a local weather model, and the measured main beam 
efficiency was used to convert the data to units of main beam temperature, Tmb- The two parts of 
the in-band frequency switched data were aligned and averaged, weighted by the inverse square of 
their individual Tsys- The data were then converted to AlPqfj SDFITS format using the GBT local 
utility idlToSdfits. In AIPS, the data were combined and gridded using the DBCON and SDGRD 
procedures. Finally, the data cubes were written to FITS files using FITTP. 



^GBTIDL is an interactive package for reduction and analysis of spectral line data taken with the GBT. 
■^The NRAO Astronomical Image Processing System 
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2.2. Australia Telescope Compact Array 

Maps of NH3 (1,1) and (2,2) inversion line emission of the Ophiuchus Bl, B2, C and F Cores 
were made over two separate observing runs at the Australia Telescope Compact Array (ATCA). 
The ATCA is located near the town of Narrabri in New South Wales, Australia, and consists of six 
antennas, each 22 m in diameter. Five of the six antennas are movable along the facility's east-west 
line and small north line, while the sixth antenna is permanently placed along the east-west line 
at a distance of 6 km from the other antennas. An 8 MHz bandwidth with 1024 channels was used, 
which provided a spectral resolution of 7.81kHz (0.1kms~^ at 23.694 GHz). This configuration 
enabled in-band frequency switching for the observations. The Oph B Cores were observed over 
four 9-hour tracks of the array (August 5 - 8, 2004), and the Oph C and F Cores were observed 
over three 9-hour tracks of the array (May 5 - 7, 2005). Both sets of observations were done with 
the array in its H168 configuration. This is a compact, hybrid configuration, where three of the 
movable antennas are placed along the east-west line and two antennas are located on the north 
spur. Baselines ranged from 61.2 m to 184.9 m (~ 4.7 kA - 14.2 kA) with five antennas. At 23.7 GHz, 
these observations provided a primary beam (field-of-view) of 2'. 

Maps were made of the cores using separate pointings spaced by &n = (2 / \/3) A / 2D ~ 1.1' 
at 23.7 GHz for Nyquist sampling on a hexagonal grid. Table [3] gives the number of individual 
pointings required to cover each Core, the multiple-beam overlap area observed, and the final rms 
sensitivity achieved towards each Core for both the ATCA and the VLA observations (described 
below) . 

Observations cycled through each individual pointing between phase calibrator observations 
to maximize uv-coverage and minimize phase errors. The phase calibrator, 1622-297, was observed 
every 20 minutes, and was also used to check pointing every hour. Flux and bandpass calibration 
observations were performed every shift on the bright continuum sources 1253-055, 1934-638 and 
1921-293. 



The ATCA data were reduced using the MIRIAD data reduction package (jSault et al.lll995l ). 
The data were first fiagged to remove target observations unbracketed by phase calibrator mea- 
surements, data with poor phase stability or anomalous amplitude measurements. The majority 
of the data were good, as the weather was stable during the observations. Much of the data from 
baselines involving the 6 km antenna, however, were fiagged due to poor phase stability. The band- 
pass, gains and phase calibrations were applied, and the data were then jointly deconvolved. First, 
the data were transformed from the spatial frequency {u, v) plane into the image plane using the 
task INVERT and natural weighting to maximize signal to noise. The data were then deconvolved 
and restored using the tasks MOSSDI and RESTOR to remove the beam response from the image. 



MOSSDI uses a Steer-Dewdney-Ito CLEAN algorithm (jSteer et al.lll984l ). The cleaning limit was 
set at twice the rms noise level in the beam overlap region for each object. Clean boxes were used 
to avoid cleaning noise in the outer regions with less beam overlap. Applying natural weighting 
provided a final synthesized beam of ~ 8" x 10" FWHM. 
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2.3. Very Large Array 

Maps of NH3 (1,1) and (2,2) emission were made at the Very Large Array (VLA) near Socorro, 
NM, USA over the period of 2007 January 20 through 2007 February 11. Nine observing shifts were 
ahotted to the project in the array's DnC configuration, each five hours in duration covering the 
LST range 14:00 - 19:00. The DnC configuration is a hybrid of the most compact D configuration 
with the next most compact C configuration. This setup ensures a more circular beam shape for 
southern sources hke Ophiuchus while retaining the sensitivity of the D configuration. 

For these observations, we used a correlator setup with two IFs, each with a bandwidth of 
3.125 MHz with 24.414kHz spectral resolution (0.3kms~^). While not providing the same high 
spectral resolution obtained at the ATCA, this setup enabled simultaneous observations of the (1,1) 
and (2,2) lines, and allowed the main hyperfine component and the two middle satellite components 
of the (1,1) line to be contained within the band. Mosaic maps were made with Nyquist-spaced 
(~ 1.0' at 23 GHz) individual pointings. Observations cycled through pointings between phase 
calibrator observations. Table [3] gives the pointing and rms sensitivity information for each Core. 

During these observations, several antennas in the array had been upgraded as part of the 
Expanded VLA (EVLA) project. Online Doppler tracking at this time was not yet available, and 
the observations were thus obtained in a fixed frequency mode, with line sky frequencies calculated 
using the NRAO's online Dopset tool. The observing frequencies were updated frequently, with 
phase calibrator observations on 1625-254 before and after a frequency change to avoid phase jump 
problems. Bandpass and absolute flux calibration were done for each shift using observations of 
1331+305. 

The data were checked, flagged and calibrated using the NRAO Astronomical Image Processing 
System (AIPS), following the procedures outlined in the AIPS Cookbooli^. In addition, special 
processing was required to account for differences in bandpass shape and antenna sensitivity between 
the VLA and EVLA antennas. System temperatures were lacking for EVLA antennas, so the VLA 
back-end Tsys values were used when reading the data into an AIPS uf-database. A bandpass table 
was then created from the line dataset and applied to a spectrally averaged 'channel 0' dataset. 
The normal VLA calibration steps were then followed. The calibrated uv flies were then written to 
FITS format and imported to MIRIAD, where the data were deconvolved and restored. Applying 
natural weighting and a taper of 8" x 6" provided a flnal synthesized beam of ~ 10.5" x 8.5" 
FWHM. 



''http://www.aips.nrao.edu/CookHTML/CookBook.html 
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2.4. Combining Single Dish and Interferometer Data Sets 

Since none of the antennas in the ATCA and the VLA can act as single radio dishes, there 
is necessarily an upper limit to the size of structure to which each is sensitive. This upper scale 
limit is dependent on the shortest spacing between two antennas in the array, and the missing 
information is thus referred to as the short- or zero-spacing problem. Mosaicing helps to recover 
short spacing information. For complex sources with emission on many scales, however, combining 
interferometer data with single-dish observations provides more complete coverage of the uv plane 
and thus creates a more accurate representation of the true source emission structure. Ideally, the 
single dish diameter should be larger than the minimum interferometer baseline to ensure maximal 
overlap in the uv-plane and determine accurately flux calibration factors between the observations. 

Data from each interferometer were first combined separately with the single dish observations. 
The GBT data were regridded to match the interferometer data in pixel scale and spectral resolu- 
tion. The data were then converted to units of Jybeam"^ for combination with the interferometer 
data using the average beam FWHM measured at the GBT during the observations. Combina- 
tion of the data was done using MIRIAD's IMMERGE task. IMMERGE combines deconvolved 
interferometer data with single-dish observations by Fourier transforming both datasets and com- 
bining them in the Fourier domain, applying tapering functions such that at small spacings (low 
spatial frequencies) the single-dish data are more highly weighted than the interferometer data, 
while conversely at high spatial frequencies the interferometer data are weighted more highly. The 
flux calibration factor between the two datasets was calculated in IMMERGE by specifying the 
overlapping spatial frequencies between the GBT and the interferometers. We took the overlap 
region in both cases to be 35m - 100m (2.7kA - 7.7kA), yielding a flux calibration factor of ~ 1.4 
between the GBT and ATCA datasets, and ~ 1.0 between the GBT and VLA datasets in the NH3 
(1,1) line emission. These factors were then applied to the NH3 (2,2) datasets. The final resolution 
of the combined data is the same as that of the interferometer data. 

To combine all three datasets, the ATCA and VLA data were imaged together using the 
INVERT task, applying natural weighting and taper as described for the VLA imaging. The 
interferometer data were then cleaned and combined with the GBT data as described above. The 
overlap region was taken to be 35 m - 100 m, yielding a flux calibration factor of ~ 1.3. By convolving 
the combined data to match the 32" resolution of the GBT data, we estimate the total flux of the 
combined image recovers nearly all 98 %) the flux in the single-dish image. These data were 
used to identify structures in the NH3 data cubes using an automated structure finding algorithm 
(sec §3). The data were also tapered to a slightly lower resolution of 15" FWHM to provide higher 
signal-to-noise ratios for a multiple component hyperfine line fitting routine. 
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Results 



3.1. Comparison with submillimeter dust continuum emission 



We first discuss the NH3 (1,1) intensity in the combined datasets, and compare the distribution 
of NH3 emission and 850 continuum emission in the Oph Cores as shown in Figure [TJ Figures 
[2^, [3^ andH^ show the combined NH3 (1,1) hue emission at ~ 10.6" x 8.5" FWHM resolution. 
The emission has been integrated over the central hyperfine components in the Oph B, C and F 
Cores, respectively, with a clip of ~ 2x the map rms noise level. (Since the outer edge of the 
combined maps have higher rms noise levels, integrating only over, i.e., the 'main component' of 
NH3 reduced the amount of signal included in the noisy outer regions.) The respective 850 ;um 
emission for each Core at ~ 15" FWHM resolution, i.e., lower than the resolution of the combined 
NH3 data, is shown in Figures [2]3, [8)3 and Bb- We also show loc ations of submillimeter clumps 
identified with the 2D version of clumpfi nd (IJorgensen et al.ll2008l ). In Figure [2)3, we additionally 
label the ~ 1 Mq continuum object MM8 (|Motte et al.lll998l ). 



Figures [21 [3] and [4] also show locations and labels of 'cold' YSOs (based on bolometric tem- 
peratures derived from fittin g their spectral ene rgy distributions) detected and classified through 
Spitzer infrared observations (lEnoch et al.ll2008l ) overlaid on the submillimeter and NH3 emission. 
The objects plotted were all identified as Class I protostars (no Class protostars have been asso- 
ciated with Oph B, C or F). Oph B2 is associated with three YSOs. Of these, two are previously 
known (IRS45/GY 273 and IRS47/GY279). Based on association with a continuum emission object 
([Enoch et al.ll2008l ) and observed infrared colours ([Jorgensen et al.[[2008l ) , either three or two proto- 
stars in Oph B2 are embedded in gas and dust. One additional Class I protostar is located between 
Bl and B2, and is only associated with diffuse NH3 emission. Oph Bl and B3 appear starless. No 
embedded pr otostars are associated with Oph C by [Enoch et al.[ ([2008[ ). but a "Candidate YSO" is 
identified by ([Jorgensen et al.l[2008[ ) 30" south of the Core continuum peak (R.A. 16:26:59.1, Dec. 
-24:35:03). Based on the different classifications by the two papers, the significant offset from the 
continuum emission peak, and the lack of any clear influence on the gas in our data, we will discuss 
Oph C assuming it is not associated with a deeply embedded object. Four protostars are associated 
with Oph F, all of which have been previously identified (see Figure [4[ for object names). Three 
are likely embedded in the Core. 

Peaks of NH3 integrated intensity can be used to surmise the presence of 'objects', but such 
identifications can ignore any differences in velocity between adjacent den se gas. We therefore 
used the 3D version of the automated structure- finding routine clumpfind ([Williams et al.[[l994l ) 
to identify distinct NH3 emission objects in the combined NH3 (1,1) data cube, which we will 
henceforth call "NH3 clumps", clumpfind uses specified brightness contour intervals to search 
through the data cube for distinct objects identified by closed contours. The size of the identified 
clumps are determined by including adjacent pixels down to an intensity threshold, or until the 
outer edges of two separate clumps meet. The clump location is defined as that of the emission peak. 
CLUMPFIND was used only on the main emission component, where multiple hyperfine components 
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are sufficiently blended to present effectively a single line given the 0.3kms~^ spectral resolution 
and typically wide line widths found. The standard interval of 2x the data rms between contours, 
with a slightly larger threshold value worked well to separate distinct emission components in all 
clumps. A lower limit of IK with intervals of 0.4 K, ~ 2x the cube rms, identified separate 
emission peaks sufficiently for Oph B and C. Oph F was sufficiently fit using a lower limit of 1.2 K 
and intervals of 0.65 K, as the rms of the combined data was slightly higher. Even so, some identified 
clumps appeared by eye to be noise spikes at the map edges, and these were culled from the final list, 
as well as any clumps which contained fewer pixels than in the synthesized beam. Table [5] lists the 
locations, FWHMs, effective radii and peak brightness temperatures for CLUMPFiND-identified cores 
in Oph B, C and F, and the clump locations are overlaid on Figures [2] - HI The centroid velocities 
and line widths of the NH3 clumps were determined through spectral line fitting (described further 
in §4), which provided more accurate measures of vlsr and Aw given the hyperfine structure of 
the NH3 lines. 

In Oph B, Figure [2] shows that the large-scale structure of the Core is similar in both line 
and continuum emission, but some significant differences are also apparent. For example, the 
integrated line emission displays a less pronounced division between the Bl and B2 Cores seen in 
the continuum, and indeed shows significant filamentary structure in the region connecting Bl and 
B2 where little to no continuum emission is observed. Similarly, strong line emission is present in the 
south-eastern edge of B2 where the continuum map shows relatively little emission. Additionally, 
while B2 is the stronger continuum emitter, Bl is significantly brighter in integrated line emission. 
Strong NH3 emission in the northern half of Bl extends beyond the continuum contours, and 
becomes significantly offset from the bulk of the continuum emission to the east . The NH3 (1,1) 
observations of the Oph B core region also reveal a narrow-lin e emission peak n orth of the western 
edge of B2, which is coincident with a DCO"*" object, Oph B3 (iLoren et al.lll99d ). B3 is just visible 
at the lowest contour in the NH3 integrated intensity map but is not visible in the continuum map. 

Furthermore, although the peaks of continuum emission and line emission are often located in 
the same vicinity in Bl and B2, the brightest continuum peaks and the integrated line emission 
maxima are typically non-coincident. Overall, the mean separation between an NH3 integrated 
intensity peak and the nearest 850 /im continuum peak in Bl and B2 is ~ 22" (2600 AU), or 
~ 2x the NH3 FWHM resolution. Fifteen CLUMPFiND-identified NH3 clumps are identified in the 
combined Oph B map, with an average minimum separation between NH3 clumps of 47" (5600 AU, 
or 40" (4800 AU), if Oph B3 is not included). The mean minimum distance between NH3 clumps 
and submillimeter continuum clumps is 44" (5300 AU), or ~ 4x the NH3 FWHM resolution. Only 
five of the fifteen NH3 clumps are located within 30" (3600 AU) of a submillimeter continuum 
clump. 

No protostars in Oph B are found at positions of NH3 (1,1) integrated intensity maxima 
nor are they associated with identified NH3 clumps. One protostar, IRS47/GY279, is located 
south of the NH3 clump we identify as B2-A6. The offset between the NH3 clump peak and the 
protostar, however, is ~ 30", or approximately 3x the angular resolution of the combined NH3 
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data. A second protostar, IRS45/GY273, is coincident with a submillimeter continuum peak but 
has httle associated NH3 emission. E ither of these pr o tostar s may be the source of an east-west 



ahgned outflow recently proposed by iKamazaki et al.l (|2003l ) from CO observations of B2. The 



third, previously unidentified protostar, seen west of NH3 clump B2-A4, is located within a narrow 
(~ 30") NH3 integrated intensity minimum between the B2 core and the filament connecting Oph 
Bl and B2. 

Figure [3] shows that Oph C has more similar overall structure when traced by the continuum 
and integrated NH3 (1,1) line emission than Oph B. Extended emission in Oph C is elongated along 
a southeast-northwest axis and contains a single integrated intensity peak. A thin (~ 30") filament 
of faint emission extends to the north of the central peak. The submillimeter continuum emission 
is largely coincident with the integrated NH3 contours, but the continuum emission peak is offset 
to the NH3 integrated intensity peak by ~ 20". clumpfind separates the central NH3 emission 
into two cores, C-Al and C-A3, and finds a third object, C-A2, at the tip of the northern emission 
extension. C-Al and C-A3 are found to the northwest and southeast (30" offset and 15" offset, 
respectively) of the centres of both the NH3 integrated intensity and of the continuum emission. 
Continuum emission also extends in the direction of the faint NH3 northern extension, but there 
is no secondary peak present. 

Like Oph C (but unlike Oph B) , Figure H] shows that Oph F also has very similar structure 
when traced by either the submillimeter continuum emission or integrated NH3 intensity. Unlike 
both Oph B and C, the CLUMPFiND-identified NH3 clumps are coincident with the integrated in- 
tensity peaks. F-A3 is nearly coincident (within a beam FWHM) with a submillimeter continuum 
clump and is additionally coincident with an embedded protostar, IRS43/GY265. F-A2 is associ- 
ated with a second embedded protostar (CRBR65) and continuum emission, but not an identified 
continuum clump. A thin filament (~ 15") extends to the northwest and the third NH3 clump, 
F-Al, which is also coincident with extended continuum emission but no identified clump. A third 
embedded protostar (IRS44/GY259) in the north-east is coincident with a continuum peak, but 
has no associated NH3 emission. A fourth protostar, in the south-west, may be coincident with 
some unresolved NH3 emission, but is located in a section of the map with larger rms values and 
consequently the small integrated intensity peak seen at that location may be simply noise. 

The discrepancies between NH3 and submillimeter continuum emission in Oph B are in con- 
trast to earlier findings of extremely hig h spatial correlation between the two gas tracers in isolated, 



low-mass starless clumps. For example. iTafalla et al.l (j2002l ) found that both the millimetre contin- 



uum and integrated NH3 (1,1) and (2,2) line intensity were compact and centrally concentrated in 
a survey of five starless clumps, including L1544 in the Taurus molecular cloud. In these clumps, 
the integrated intensity maxima of both transitions are approximately coincident (within the 40" 
angular resolution of the NH3 observations) with the continuum e mission peaks. This same coinci- 



dence between NH3 and millimeter continuum was found in B68 (iLai et al.ll2003l ). When observed 
at higher angular resolution, NH3 emission in LI 544 remained coincident with the continuum but 
the line integrated intensity peak was offset by ~ 20". This offset was explained, however, as being 
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due to the the NH3 emission becoming optically thick (jCrapsi et al.ll2007l ). 



It is also possible th at different methods of id entifying structure in molecula .r gas (such as the 



GAUSS CLUMP method of IStutzki &: GuestenI ll990|, or using dendrograms as in iRosolowskv et al 



2008bl . for example) would create a different 'core' list than presented here. We have additionally 
compared our result s in Oph B with the locations of millimeter objects identified using multi- 
wavelet analysis by iMotte et alj (|l998l ). and find a smaller yet still significant mean minimum 
distance of 27" between NH3 clumps and millimeter objects. Given the severe positional offsets in 
some locations between the NH3 emission and submillimeter continuum, it is unlikely that different 
structure- finding methods would provide substantially different results. 



3.2. Single Dish C2S and HC5N Detections 

The GBT observations of C2S 2\ — Ig emission only resulted in single, localized detections 
in Oph Bl and Oph C, while only Oph C had a single, localized detection in HC5N (9-8). The 
C2S emission in Bl was confined to a single peak at its southern tip. The single-dish spectra of 
all observed species at the C2S peak locations in Bl and C are presented in Figure [5] (note that 
in Figure [5l the NH3 line is so narrow in Oph C that we are detecting the hyperfine structure 
of the (1,1) line). The integrated intensity GBT maps of all molecules observed in Oph C are 
shown in Figure [6l Within the ~ 30" resolution limits of the GBT data, the NH3, C2S and HC5N 
spectral line integrated intensity peaks overlap with the local 850 //m continuum emission peak in 
Oph C. We fit the spectra of the two C2S and single HC5N detections with single Gaussians to 
determine their respective vlsRj line width Av, and peak intensity of the line in Tmb units. We 
additionally fit a 2D Gaussian to the integrated intensity maps to determine the FWHMs of the 
emitting regions. We find a beam-deconvolved FWHM = 9200 AU and 5800 AU for the C2S and 
HC5N emission, respectively, in Oph C. The C2S emission in Oph Bl is elongated, with a FWHM 
= 9400 AU in R.A. but only 3700 AU in Dec. for an effective FWHM = 5100 AU. The results of 
the Gaussian fitting are listed in Tabled! See §4.4.3 for further analysis of C2S and HC5N. 



4. NH3 Line Analysis 

4.1. NH3 Hyperfine Structure Fitting 

The metastable J = K rotational states of the symmetric-top NH3 molecule are split into 
inversion doublets due to the ability of the N-atom to quantum tunnel through the hydrogen atom 
plane. Quadrupole and nuclear hyperfine effects further split these inversion transitions, resulting in 
hyperfine structure of the (J, K) = (1,1) transition, for example, containing 18 separate components. 
This hyperfine structure allows the direct determination of the optical depth of the line through the 
relative peaks of the components. Additionally, since transitions between A'-ladders are forbidden 



- 12 - 



radiatively, the rotational temperature describing the relative populations of two rotational states, 
such as the (1,1.) ^.nd (2,2) transitions, can be used to determine directly the kinetic gas temperature 



(|Ho fc Towneslll983l ). 



For a given (J,K) inversion transition of NH3 in local thermodynamic equilibrium (LTE), the 
observed brightness temperature as a function of frequency v can be written as 

TX,{J,K) = m4B ^ {J{TMK)) - J{ng)){l - eM-ru{J,K))) (1) 



L ex,v 



assuming the excitation conditions of all hyperfine components are equal and constant (i.e., Te: 
Tex)- Here, tjmb is the main beam efficiency, <I> is the beam filling factor of the emitting source, 
Tex is the line excitation temperature, T^g = 2.73 K is the temperature of the cosmic microwave 
background, and J{T) = (hv / k)[exp{hi^ / kT) — 1]~^. The line opacity as a function of frequency, 
Tu, is given by 



= To aj exp (^-A In 2 ^ 



(2) 



where is the total number of hyperfine components of the (J,K) transition (N = 18 for the 
(1,1) transition and = 21 for the (2,2) transition). For a given j'*'' hyperfine component, aj is 
the emitted line fraction and Vj is the expected emission frequency. The observed frequency of 
the brightest l i ne com ponent is given by fQ, with a FWHM Au. Values of aj and Vj were taken 



from iKukolichI (jl967l ). Here, we assume $ = 1. If the observed emission does not entirely fill the 
beam, the determined T^x will be a lower limit. In regions where the emission is very optically thin 
(r << 1) there is a degeneracy between r and T^x and solving for the parameters independently 
becomes impossible. We restricted our analysis to regions where the NH3 (1,1) intensity in the 
central component is greater than 2 K, which corresponds roughly in our data to a signal-to-noise 
ratio of 8 — 10 in the main component and ~ 2 — 3 in the satellite components. With this restriction, 
we also find r > 0.5 throughout the regions discussed. 

To improve the signal-to-noise ratio of the data and to match the resolution of the 850 ^m 
continuum data, we first convolved the combined data to a final FWHM of 15" (from 10.6" x 
8.5"), and then binned the convolved data to 15" x 15" pixels. Assuming Gaussian profiles, the 
18 components of the NH3 (1,1) emission line were fit simultaneously using a chi-square reduction 
routine custom written in IDL. The returned fits provide estimates of the line centroid velocity 
(visr), the observed line FWHM (Avobg), the opacity of the line summed over the 18 components 
(r), and [J (Tex) — J{Tf,g)]. The satellite components of the (2,2) line are not visible above the 
rms noise of our data. These data were consequently fit (again in idl) with a single Gaussian 
component. 



The line widths determined by the hyperfine structure fitting routine are artificially broadened 
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by the velocity resolution (0.3 km s ^) of the observations. To remove this effect, we subtract in 
quadrature the resolution width, Au^es; from the observed line width, Avobsi such that Avune = 
Y''Az^^^^"^-Au^. In the following, we simply use Av = Avune for clarity. The limitations of the 
moderately poor velocity resolution are discussed further in Appendix [Bl but do not significantly 
impact our analysis. In regions where lines are intrinsically narrow, such as Oph C and parts of 
Oph F, the derived line widths may be overestimated by up to ~ 20 — 35 %. 

The uncertainties reported in the returned parameters are those determined by the fitting rou- 
tine, and do not take the calibration uncertainty of ~ 8% into account. The calibration uncertainty 
affects neither the derived parameters that are dependent on ratios of line intensities, such as the 
opacity and kinetic temperature Tk, nor the uncertainties returned for vlsr or Av. The excitation 
temperature, however, as well as the column densities and fractional NH3 abundances discussed 
below (see §4.4) are dependent on the amplitude of the line emission, and are thus affected by the 
absolute calibration uncertainty. 

Table [6] lists the mean, rms, minimum and maximum values of Visr, Av, r and T^x found in 
each of the Cores using the above restrictions for the combined NH3 line emission. In the following 
sections, we describe in detail the results of the line fitting and examine the resulting line centroid 
velocities and widths, as well as T^x and r. In addition, we use the fit parameters to calculate the 
gas kinetic temperature (Tk), non-thermal line widths {(Jnt), NH3 column density (iV(NH3)), gas 
density (n(H2)) and NII3 abundance (X(NH3)) across all the cores, as described further in §4.2, 4.3 
and 4.4. Table [7] summarizes the mean, rms, minimum and maximum values of Tk, ctnt, -^(NH3), 
n(Il2) and X(NIl3) for each Core. For each NII3 clump. Table [8] summarizes the mean, rms, 
minimum and maximum values of all determined parameters (means were obtained by uniformly 
weighting each pixel). 



4.2. Line Centroids and Widths 

Figures [7^1, [8^, and [9^ show maps of vlsr of the fitted NII3 (1,1) line in Oph B, C and F 
respectively from the combined, smoothed and regridded data. These maps reveal that although 
variations of vlsr are seen within the Cores, they are not that kinematically distinct from each 
other. For example, only 0.35 km s~^ (< the mean Av) separates the average line-of-sight velocity 
in Oph B from Oph F . This result agrees with the ID velocity dispersion of ~ 0.36 km s~^ found 



by I Andre et al.l (|2007l ) through N2H'^ observations of the Oph cores. 



In Oph B (see Figure [7^), the vlsr of NII3 emission has little internal variation, with a mean 
Vlsr = 3.96 kms~^ and an rms of only 0.24kms~^. An overall gradient is seen across Oph Bl 
and B2, with smaller vlsr values (3.2 kms^^) at the southwest edge of Bl increasing to 4.6kms^^ 
at the most eastern part of B2. Correspondingly, Bl has a characteristic velocity somewhat less 
than the average (3.8kms~^ itO.lkms"^), while the B2 vlsr is slightly greater (4.1 it 0.2 kms~^). 
The filament connecting Bl and B2 is kinematically more similar to B2, but there is no visible 
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discontinuity in line-of-sight velocity of the lines. B3 has the lowest vlsr in Oph B, with an 
average velocity of 3.3kms~^ it 0.2kms~^, or 0.7kms~^ less than the average of the group. This 
large difference, greater than that between the mean vlsr values of Oph B and F, suggests that B3 
may not be at the same physical distance as the rest of Oph B. The change in velocity between B2 
and B3 occurs over a small projected distance (~ 30", or ~ 3600 AU). There is some indication that 
B2 and B3 may overlap along the line of sight, as the determined vlsr values in B2 immediately 
south of B3 are less than those to the east and west, as might be expected if lower velocity emission 
is also contributing to the line at that location. (As discussed further below, the A?; line widths in 
this area are larger than the average, as would be expected if unresolved emission from two different 
velocities is contributing to the observed line.) 

Oph C (see Figure EH) has a mean vlsr = 4.01 kms^^ with an rms of only 0.07 kms~^. A 
small velocity gradient of ~ 0.4kms~^ is evident in Oph C, with a minimum line-of-sight velocity 
of ~ 3.7kms~^ in the southeast, increasing to a maxi mum of ^ 4. 1 kms~ ^ in the northwest. A 



similar gradient was noted in N2H^(l-0) observations by I Andre et al.l (120071 ). and may be indicative 
of rotation. The identification of two NH3 clumps in the region, however, could also be indicative 
of two objects with slightly different vlsr- A slight decrease in vlsr values is seen in the NH3 
extension to the north, with C-N2 associated with emission at a slightly lower vlsr than the mean. 

Oph F (see Figure [9^) has a mean vlsr = 4.34 km s~^ and an rms of only 0.26 km s~^ over 
the region containing the bulk of the NH3 (1,1) emission and two of the three identified NH3 
peaks, despite the presence of four protostars. The filament extending towards the northwestern 
NH3 peak gradually increases in vlsr^ but only by ~ 0.2kms~^. With higher sensitivity, single- 
dish data show that outside this area vlsr drops to ~ 3.7 km s^^. There is clear evidence for 
two velocity components alo ng the line o f sight at the peak position of F-Al, with a secondary 



component at ~ 3.7kms~ . lAndre et al.l (|2007l ) also find two velocity components near F-Al in 
N2H+observations. For the brighter component, their N2H+ data agree with our NH3 data, but for 
the secondary component they find a higher vlsr = 4.1kms~^. Some blue asymmetry is found in 
the line profiles of F-A2 and F-A3, but this is more likely due to the complicated velocity structure 
of the core rather than infall motions (comparison with an optically thin tracer at this position is 
necessary to confirm infall). 

In summary, vlsr varies little across any of the Cores (rms < 0.24kms^^), and additionally 
varies little between them (maximum mean difference ~ 0.38 km s^^). Some small gradients were 
found in the larger Cores (i.e., on scales larger than the individual NH3 clumps) which may be 
indicative of rotation. 

Figures Eh, Eb, and [9]o show the line Af for Oph B, C and F obtained from the combined, 
smoothed and regridded data. As stated above, the line widths have been corrected for the reso- 
lution of the spectrometer (0.3kms~^). In the few cases where the returned FWHM from the fits 
is similar or equal to the resolution, we set the corrected FWHM to the thermal line width (see 
Appendix A for further discussion). 
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Line widths range from Af < O.lkms" to ~ 1.2kms~ in both Oph B and F, with a 
minimum hne width of 0.08 km in B3 and a maximum of 1.37 km s^^ in B2. Line widths in 
Oph C range from 0.11 kms~^ to 0.70 kms~^. The rms variations Av and vlsr in Oph B and C 
are similar (0.2kms~^ and O.lkms"^, respectively), while in Oph F the rms in Av is much larger 
than the rms in vlsr (i-e., 0.3kms~^ compared to 0.12 kms^^, respectively). 

The extended emission in Oph B is dominated by highly non-thermal motions (mean Av = 
0.83 kms~^), shown in Figure[7)3. Several localized pockets of narrow line width are found embedded 
within the more turbulent gas. The single Oph B3 clump (B3-A1) and B2-A7 are both characterized 
by extremely narrow Av (i.e., 0.08 kms~^ and 0.33 kms~^ respectively). As mentioned above, the 
line emission broadens from B3 to B2 over only ~ 30" (3600 AU) to At; ~ 1.4kms-^ possibly 
due to line blending along the line of sight if B2 and B3 overlap in projection at these locations. 
We find small Av ~ 0.5kms~^ towards the south-eastern edge of the mapped region in Bl near 
NH3 clumps B1-A3 and B1-A4, but the line width minimum does not coincide with either core. 
Another region of low Av ~ 0.6kms~^ is coincident with the NH3 clumps B2-A1 and B2-A2. A 
final At; minimum, Af ~ 0.4kms~^ is found between the two eastern protostars in B2. Note that 
the protostars in Oph B are all associated with smaller Af than the average value for the core, but 
none are coincident with a clear local minimum in line width. 

The maximum line width in Oph C, Af = 0.70 km s~^, is half that found in Oph B. Figure 
[Hb shows that most of the emission in C is narrow, with a minimum Af = 0.11 km s~^, similar 
to the extremely narrow lines found in B3. These narrowest line widths are found centered on 
NH3 clump C-A2 in a band perpendicular to the elongated direction of Oph C and are coincident 
with the highest velocity emission, and are thus not coincident with the NH3 integrated intensity 
maximum nor the continuum emission peak. Curiously, if C-A2 and C-A3 are indeed physically 
distinct clumps, we would expect the broadest lines between them due to overlap, but instead find 
narrow lines at this location. The largest line widths are found at the edges of the integrated 
intensity contours. 

Oph F is characterized by moderately wide line emission more similar to that found in Oph B, 
with a mean Af = 0.63 kms~^ (see FigurelHb). Line widths associated with the central NH3 clumps 
F-A2 and F-A3 are smaller than the mean (i.e., 0.6kms~^ and 0.3kms~^ respectively), while clear 
minima in line width (Af ~ 0.33 kms~^) are found at the locations of the two central protostars. 
Line widths along the Core extending to the northwest integrated intensity peak broaden to larger 
values (~ 1.2kms~^), but at the tip NH3 clump F-Al is associated with Af = 0.4kms~"'^ in a 
single 15" pixel. 

Overall, the observed NH3 Af in the Cores are generally large, excepting Oph C. We also 
find regions of localized narrow line emission, some of which are associated with NH3 clumps. 
Line widths near protostars tend to be smaller than the mean values, but only in Oph F are the 
protostars coincident with clear minima in Af . 
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4.3. Kinetic Temperatures and Non-Thermal Line widths 



Following iMangum et alj (|l992l ). among others, we use the returned r, A?; and line brightnesses 
of the NH3 (1,1) and (2,2) lines in each pixel to calculate the kinetic temperature of the gas. Details 
of our calculations can be found in Appendix [Al Propagating uncertainties from our hyperfine 
structure fitting routine gives typical uncertainties in Tk of ~ 1 K. The mean, rms, minimum and 
maximum kinetic temperatures are given for Oph B, C and F in Table [7l Figures [7b, Eb and [9b 
show the kinetic temperatures calculated across the Oph Cores. 

In Oph B, we find a mean Tk = 15.1 K with an rms variation across the entire core of only 
1.8 K. Most NH3 emission peaks are associated with lower than average gas temperatures, but only 
a few are coincident with clear Tk minima. The lowest temperatures in the Core, Tk = 13.7 K and 
12.4 K, are found in southern Bl towards B1-A3 and B1-A4, respectively. These low temperatures 
are found at the same location as the detection of single dish C2S emission. Gas temperatures appear 
colder [Tk ~ 13 K) towards the centre of Oph B2, but the minimum temperature, Tk = 12.4 K, is 
not coincident with an NH3 clump. Instead, the lowest temperatures are found directly between 
the NH3 clumps B2-A4, B2-A5 and B2-A6, and closer to the central submillimeter clump. 

Oph C is the coldest of the observed cores, with a mean Tk = 12.8 K and a similar rms variation 
(1.6 K) as in Oph B. The central region is effectively at a single low temperature Tk = 10.6 K, with 
the lowest values found near the emission peaks C-Al and C-A3 in the northwest and southeast, 
while the gas temperature of the northern core, C-A2, is consistent with the average. 

Oph F is characterized by the highest temperatures of the observed Cores, with a mean Tk = 
16.6 K, slightly warmer than Oph B, and with a large rms variation of 3.2 K. No clear minima in gas 
temperature are observed near any of the NII3 clumps, protostars or continuum peaks identified 
in the region. The protostar associated with NII3 clump F-A2 is coincident with a temperature 
maximum in a single 15" pixel. 

The gas temperatures traced by NH3 emission in the Oph B and F Cores are consistently higher 
than those found in isolated dense clumps. For example, all five of the starless clumps surveyed by 
Tafalla et al.l (|2002l ) were found to have a constant gas temperature Tk = 10 K determined through 
an analysis similar to that done here. Two recent studies of NH3 emission in dense clumps in the 
Perseus molecular cloud and the less active Pipe Nebula at 32" resolutio n also found slightly lowe r 
temperatures than those found here, with a median Tk = 11 K in Pe rseus (iRosolowsky et al.ll2008al ) 
and a mean Tk = 13 K ±3 K for < 1 Mq clumps in the Pipe Nebula ( Rathborne et al.ll2008l . we note 
that one object in the Pipe is warmer than the typical Tk found i n Oph B and F). In a sample of 
NII3 observations towards Galactic high mass star forming regions, IWu et al.l (|2006l ) found a mean 
Tk = 19 K. The mean kinetic temperatures found in Oph B and F are also slightly gr eater than 
the median Tk = 14.7 K found in a survey of NII3 observations by iJijina et al.l (|l999l ). although 
their analysis showed that the median temperature of dense gas in clusters was significantly higher, 
Tk = 20.5 K, than in non-clustered environments where the median Tk = 12.4 K. Since L1688 is 
a clustered star forming environment, it is not unreasonable to expect temperatures higher than 
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those in more isolated regions, given the lJiiina et al.l results. 



Evidence for an extremely cold temperatm'e of 6 K [obtained through observations of H2D^" (lio~ 
111), which li kely probes dense r gas than NH3 (1,1) and (2,2), was recently found for the nearby 
Oph D Core (iHarju et al.ll2008l ). In addition, while the mean gas temperature Tk = 12K in Oph 
C, temperatures in the highest column density gas drop to 10 K, similar to temperatures found in 
the studies of isolated clumps described above. 

We find little difference in average gas temperatures between NH3 clumps and NH3 emission 
associated with submillimeter continuum emission peaks, and a mean Tk increase of only ~ IK, 
i.e., similar to our uncertainty in Tk, in gas temperatures near protostars ( but note on l y five 
protostars are associated with emission with sufficient S/N to fit the HFS). The lJijina et al.l (119991 ) 
survey found that NH3 cores not associated with identified IRAS sources (presumably protostellar 
objects) were slightly colder than those with coincident IRAS detections (12.4 K compared with 
15.0 K), but the effect was much smaller than temperature differences seen due to association with 
a cluster. 



Given the determined gas temperature Tk, we calculate the expected one-dimensional thermal 
velocity dispersion ax of the gas across the cores: 



f^T = W (3) 

Here, ks is the Boltzmann constant, ^nh3 = 17.03 is the molecular weight of NH3 in atomic 
units, and mn is the mass of the hydrogen atom. Similarly, the thermal sound speed Cs of the gas 
can be calculated using a mean molecular weight /x = 2.33. 

The non-thermal velocity dispersion a^T is given by 



'^NT = ^C^L -4 (4) 

where aobs = At)/(2V21n2). The mean, rms, minimum and maximum values for both ai\}T and the 
non-thermal to thermal velocity dispersion ratio of the gas, given by gmt / Cs, are given for each of 
Oph B, C and F in Table El 

Figures Eli, Ell and Ell show the resulting non-thermal to thermal velocity dispersion ratio over 
the cores. The mean umt I Cs values show supersonic velocities are present. At the limits of the 
velocity resolution of our data, the smallest observed line widths are consistent with motions being 
purely thermal in nature. 

We find that the mean o^t = 0.35 km s^"*^ and cJiVT / Cs = 1.5 in Oph B. Across the Core, 
we find a moderate ittvt / Cs rms of 0.4. The majority of the gas traced by NH3 in Oph B is 
thus dominated by non-thermal, mildly supersonic motions. Several NII3 clumps are associated 
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with smaller, but still transsonic, non-thermal motions. Thermal motions dominate the observed 
line widths in only two well-defined locations in Oph B which additionally coincide with NH3 
clumps: B2-A7, with gnt/cs = 0.5, and B3-A1, where the observed line width is consistent with 
purely thermal motions. Otherwise, little difference is seen between the non-thermal line widths of 
individual cores and the surrounding gas, with a mean ctkit / = 1.4kms~^ for the NH3 clumps. 

In contrast, Oph C has a mean umt = 0.14kms^^. Consequently, the mean a^T / Cs = 0.6 with 
an rms of only 0.2, showing that thermal motions dominate the observed line widths over much of 
the Core. The minimum non-thermal line width is found associated with C-Al, which is consistent 
with purely thermal motions within our velocity resolution limits. On average, non-thermal motions 
in Oph F are also similar to the expected thermal values, with a mean a^T / Cs = 0.9, but with 
a larger spread around the mean (0.6 rms) and a maximum value {umt / Cg = 2.2) similar to that 
found in Oph B {ctnt / Cs = 2.5). The lowest values are found towards F-A2 and F-A3 and the 
nearby protostars. 

The non-thermal NH3 line widths we measure are similar to those recently found for N2H+(1- 
0) emission in the Cores at larger angular resolut ion 26"), wher e the mean gmt/cs = 1-6 it 0.3 



in Oph B, 0.9 ±0.2 in Oph C, 1.5 ±0.8 in Oph F (lAndre et al.ll2007l l. a.nd are less than t hose found 
in DCO+ emission {unt/cs ~ 2 in B2, ~ l - 1.5 in Bl, B3, C and F (|Loren et al.lll99nl ). 



4.4. Column Density and Fractional Abundance 

44.1. NH3 

Given Av, r and Tg^ from the NH3 (1,1) line fitting results, we calculate the column density 
of the upper level of the NH3 (1,1) inversion transition. We then calculate t he NH3 partition 



functio n (given Tk) to determine the total column density of NH3 following iRosolowskv et al 



(j2008al ). Relevant equations are given in Appendix lAl 



We also calculate the H2 column density, A^(H2), per pixel in the Cores from 850 fim continuum 
data using 



iV(H2) = Su/[0.mfJ-mHK^B^{Td)], 



(5) 



where Siy is the 850 //m flux density, Qm is the main-beam solid angle, /i = 2.33 is the mean 
molecular weight, ttt-h is the mass of hydrogen, Ki, is the dust opacity per unit mass at 850 fim, and 
By iTii) is the Planck fu nction at the dust temperat ure, Trf. We take Ky = 0.018 cm^ g~^, follow- 



ing 



Shirlev et al.l (|200d ). using the dust model from lOssenkopf &: Henninel (|1994| ) which describes 
grains that have coagulated for 10^ years at a density of 10^ cm~^ with accreted ice mantles and 
incorporating a gas-to-dust mass ratio of 100. The 15" resolution continuum data were regridded 
to 15" pixels to match the combined NH3 observations. The dust temperature per pixel was 
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assumed to be equal to the gas temperature derived from HFS line fitting of the combined 
NH3 observations. This assumption is expected to be good at the densities probed by our NH3 
data ( n > 10^cm~^), when therm al coupling between the gas and dust by collisions is expected to 



begin ([Goldsmith &: Langeiill978l ). If these temperatures are systematically high, however, then 



the derived N(Ii2) values are systematically low. There is a ^ 20% uncertainty in the continuum 



flux values, and estimates of can additionally vary by ~ 3 (IShirley et al.ll200Cl ). Our derived H2 



column densities consequently have uncertainties of factors of a few. 

Due to the chopping technique used in the submillimeter observations to remove the bright 
submillimeter sky, any large scale cloud emission is necessarily removed. As a result, the image 
reconstruction techniqu e produces negative features around strong emission sources, such as Oph B 



(j Johnstone et al.ll2000al ). While the flux density measurements of bright sources are likely accurate, 
emission at the core edges underestimates the true column. We thus limit our analysis to pixels 
where > O.lJybeam"^, though the rms noise level of the continuum map is ~ 0.03 Jy beam~^. 
For a dust temperature = 15 K, this flux level corresponds to A^(H2) ~ 6 x 10^^ cm^^. 

Using the calculated H2 and NH3 column densities, we have calculated per pixel the fractional 
abundance of NH3 relative to H2, X(NH3) = iV(NH3) /iV(H2) for each Core. The results of these 
calculations are shown in Figures [71 [8| and [U which show the H2 column density derived from 
submillimeter continuum data, A^(H2) and the fractional NH3 abundance, X(NH3). The mean, 
rms, minimum and maximum of the derived column density and fractional abundance in each Core 
are given in Table [71 while specific values for identified NH3 clumps are listed in Table [8l The 
A^(H2) and consequently the X(NH3) uncertainties given in Table [8] include the ~ 20% uncertainty 
in the submillimeter continuum flux values only; uncertainty in Ki, is not taken into account. 

Oph B has a mean NH3 column density of 2.2 x lO^^cm"^, with the highest A^(NH3) values 
(maximum iV(NH3) = 4.8 x 10^^ cm~^) found in Bl. Two peaks in NH3 column density are found 
in Bl which correspond closely with the integrated intensity maxima, but are offset from the Bl 
NH3 clumps. In B2, the column density also generally follows the integrated intensity contours, 
with lower values overall than in Bl. The highest NH3 column in B2 of A^(NH3) = 4 x lO^'^cm"^ 
is found towards B2-A5. The highest opacity in B2, r = 4.7, was found associated with B2-A7, 
but the NII3 column density at this location is similar to the core average. Small column density 
increases are seen at other NH3 intensity peak locations. The NII3 extension connecting Bl and 
B2 is characterized by similar NH3 column densities to those at the edges of the Core with no 
obvious A^(NH3) maxima. 

The discrepancy between the bright NII3 and faint submillimeter continuum emission in Bl 
indicates a high NH3 fractional abundance relative to B2, with fractional abundances a factor of 
> 2 — 3 higher than the typical values of X(NH3) ~ 10~^ within B2. This is shown in Figure [7|F. 
Abundance minima are seen in B2, most notably towards B2-A6 and nearby protostars. Despite 
the higher NH3 column densities in Bl and B2, the lack of submillimeter emission in the NII3 
extension connecting the two regions suggests this connecting material has a higher fractional 
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NH3 abundance. Prominent negative features in the continuum data in this extension preclude a 
quantitative X(NH3) estimate. 

The mean and maximum NH3 column densities in Oph C are similar to those found in Oph 
B ((iV(NH3)) = 2.4 X lO^^cm-2 and iV(NH3) = 5.5 x lO^^cm-^, respectively). Oph C contains 
two A^(NH3) maxima. One is coincident with C-A3 and the second is offset to the west by ~ 15" 
from C-Al. C-A3 is correspondingly associated with a maximum in NH3 fractional abundance 
(X(NH3) = 8.8 X 10~^), but C-Al is coincident with an elongated minimum X(NH3) ~ 3 x 10~^ 
that extends along the same axis perpendicular to the long axis of the core where the smallest line 
widths were found. The highest NH3 abundances A'(NH3) ~ 12 x 10~^, are found in the northern 
extension. The mean abundance in Oph C, X(NH3) = 8.2 x 10~^, is slightly more than half the 
Oph B average. 

In Oph F, the mean A^(NH3) = 1.4 x lO^^cm"^ is less than that found in Oph B and C by a 
factor of ~ 2. The maximum A^(NH3) = 2.5 x lO^'^cm"'^ is also significantly less than the maxima 
in either B or C, and is found in the emission extending to the northwest from the two central NH3 
clumps and protostars. The fractional NH3 abundances are also low compared with B and C, with 
a mean X(NH3) = 5.4 x 10~^ and a maximum X(NH3) = 1.0 x 10^^ found near but not coincident 
with F-A3. 

Studies of isolated starless clumps have determined a wide range of NH3 abundance values for 
these objects. While in some cases different methods have been used to determine H2 column density 
values than that performed here, 'typical' observed fractional abundance values i r i cold, dense 
regions tend to be on the order of a few x lO~^ to a few xlO~ ^ ([Tafalla et al.ll2006l : ICrapsi et al 



20071 : lOhishi et al.lll992l : iLarsson et al.ll2003l : iHotzel et al.ll200ll). Abundances a s low as X(NH3) = 
7 X 10~^° and 8.5 x 10~^° have been proposed for B68 (jPi Francesco et al.l I2OO2I ) and Oph A 



(jLiseau et al.ll2003l ). respectively. The values found here agree well with previous studies. The 
wide variations of X(NH3) in the same general environment suggests dramatic differences in the 
chemical states of the Cores in L1688 (see §5). 



4.4.2. C2S and HC5N 

Similarly, we can calculate the abundance of C2S and HC5N from respective emission detected 
in the single-dish data, where 



N 



SirkuQ gl 1 
hc^ g2A~ 



ul 



27ra^[J{Tex) - J{Tbg)]T, 



ul 



(6) 



is the column density of the upper state of the observed t ransi tion (jRosolowskv et al.ll2008al ). The 
values for A^u gi and §2 were taken from iPickett et al.l (|l998l ) for each transition. Assuming the 
transitions are optically thin, the observed temperature of the line Tmb = [J(Tex) — J{Tbg)]Tui. 
The partition function Z = Yli di^^Pi^T^) then used to calculate the total column density 
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of each species as for NH3, with gi and Ei values taken from iPickett et alj (|l998l ). The column 
densities thus derived are given in Table HI The molecular column densities de rived for C2S in Bl - 
and C (A^(C2S) ^ 10^^~^^ cm~^) are sim ilar to results in young starless cores (jTafalla et al.1 120061 : 
Rosolowskv et al.ll2008al : iLai et al.ll2003l'). The iV(HC5N) results agree with p revious measurements 



i n the Taurus molecula r cloud (jCodella et alj 119971 : iBenson Mverd Il983l ) and the Pipe Nebula 
(jRathborne et alJ 120081 ) . We calculate molecular abundances as above and find X{C2S) = 3.1 x 
10^1° and 1.5 x 10"^° at the CoS emission peaks in Oph Bl and C, respectively. We further find 
an abundance X(HC5N) = 4.6 x 10~^^ at the HC5N emission peak in Oph C. 



4.5. H2 Density 



Given the determined excitation and kinetic temperatures, Tex and Tk, and assuming the 
metastable states can be approximated as a two level system, we have calculated the gas density 
n(H2) from the NH3 (1,1) transition following iHo TownesI (|l983l ). We list the mean, rms variation 
and range of densities found for each core in Table [71 Note that this density is effectively a mean 
density along the line of sight. In general, we find n(H2) ~ a few x 10^ cm^^ in all three Cores, with 
only moderate variation and no clear spatial correspondence with NH3 or continuum clumps. The 
largest n(H2) values (n(H2) ~ 8 x 10^ cm~^) were found in Oph B2 towards the central continuum 
clump MM8 (labelled in Figure [2)3). While these values agree with n(H2) estimates based on NH3 
emission in other regions, they are an order of magnitude lower than estimates of Ophi uchus clump 



densities derived frorn dust c ontinuum emission studies at similar spatial resolutions (|Motte et al 
19981 : 1 Johnstone eTallboOObl ). 



5. Discussion 

5.1. Discussion of small-scale features 

5.1.1. Correlation between NH3 clumps, NH3 integrated intensity and dust clumps 

In §3, we used clumpfind to identify objects in NII3 emission within the Oph Cores in position 
and velocity space. The returned NH3 clump locations are generally found at locations of peak 
integrated NII3 intensity, with the exception of Oph C, in which we found two distinct NII3 clumps 
offset from the NH3 integrated intensity maximum. 

In Oph B, we find poor correlation between maxima of NH3 integrated intensity and thermal 
dust continuum emission. Since most NII3 clumps are located at integrated intensity maxima, 
we hence find NII3 clumps identified through clumpfind do not correlate well with dust clumps. 
Continuum dust emission is a commonly used surrogate tracer of gas column density. The observed 
flux is a function of the dust emissivity {k^) and temperature (T^). If the dust in Bl, for example, 
was colder than that in B2, the same column of dust would produce less emission. In §4.3, we 
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calculated H2 column densities assuming the dust and gas are thermally coupled. For Oph B, 
we found the H2 column density closely followed the observed continuum emission under this 
assumption (see Figure [7^ vs. Figure [2)3). The dust and gas, however, may not have the same 
temperature. If the dust is colder than the gas by a small amount (T^ = lOK compared with 
Tx = 15K, for example), the true column density of H2 could be larger by a factor of ~ 2 along 
that line of sight. If these temperature differences occur on small enough scales, e.g., at the NH3 
clump positions, they could explain the discrepancy between the locations of dust clumps and 
NII3 clumps. Thermal coupling of gas and dust is most likely to occur, however, in the coldest 
and densest clumps, i.e., exactly where we do not find correspondence between the dust and gas 
tracers. 

A more likely cause of the offset between dust and NH3 emission is fractional abundance 
variation of NII3 in the Oph Cores. If the column densities determined from the dust emission are 
accurate, then most dust clumps are associated with X(NIl3) minima. Within Bl and B2, we find 
variations in X(NIl3) of > 2 on length scales similar to the NII3 clump sizes. Models of nitrogen 
chemistry in dense regions suggest that a long timescale, greater than the free-fall time, is required 
for molecules such as NH3 to achieve steady state val ues, but that during gravitational coll apse 
X(NIl3) begins to decrease at densities n > 10^ cm~^ ( Aikawa et al.ll2005l : iFlower et al.l 120061 ). 



The C2 S molecule is easily depleted in cold , dense environments with an estimated lifetime of a 



few X 10^ yr ( de Greg orio-Monsalvo et al 



(^Suzuki et al.i,1992 : .Lai fc Crutchei 



2000 



20061) . It is thus a go od tracer of young, undepleted cores 



Tafalla et al.ll2004l ). The detection of C2S in southern 



Oph Bl and Oph C therefore suggests that these specific locations are chemically, and hence 
dynamically, younger compared with other regions. The C2S emission detected in both Bl and 
C is coincident with or only slightly offset from the integrated NH3 intensity peak, suggesting 
significant depletion has not yet occurred at those particular locations. This conclusion is further 
bolstered by the fact that we find higher gas densities (see §4.5) in B2 than in Bl or C, and both B2 
and F are associated with embedded protostars and are therefore likely more dynamically evolved, 
i.e., denser. The higher levels of non-thermal motions found in Oph B are at odds with what is 
expected for an evolved, star forming core, however, and we discuss this further below. 



5.1.2. Comparison of NH3 clumps, submillimeter clumps and protostars 

We next compare the mean properties of the dense gas associated with the locations of NII3 
clumps, submillimeter clumps and protostars. We note that given the poor correlation between the 
NH3 clumps and submillimeter clumps and protostars, the derived physical properties (e.g., Tk 
and (Tatt) at the submillimeter clump and protostellar locations may be associated with larger-scale 
gas emission along the line-of-sight rather than the dense clump gas. 

In general, we find only small differences between the mean properties of the dense gas at 
the peak locations of the NII3 clumps, submillimeter clumps and embedded protostars. The mean 
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kinetic temperatures for NH3 clumps and submillimeter clumps are nearly equal (~ 14 K), and only 
~ 1.5 K less than the values associated with embedded protostars. This difference is not significant 
given that the uncertainties in Tk are on the order of 1 K. Conversely, excitation temperatures 
associated with embedded protostars are ~ 1.5 K lower than that of NH3 clumps and submillimeter 
clumps where T^^^ ~ 10 K, with uncertainties in T^^ also ~ IK. The line widths of submillimeter 
clumps tend to be larger than those associated with NH3 clumps by only ~ 25%. 

Some differences in mean properties between objects are notable. For example, protostars have 
associated Av and ctatt / a factor of 2 narrower than both submillimeter and NH3 clumps. Also, 
submillimeter clumps and protostars have lower fractional abundances than seen for NH3 clumps 
by a factor of ~ 2. Note, however, that only five protostars are found with NH3 emission strong 
enough to analyze, as we described in §3. For this reason, our comparison sample is limited to 
protostars that are still associated with significant amounts of gas, where conditions are likely more 
similar to those found in submillimeter clumps than for more evolved protostars. 



5.1.3. (TNT I Cs Individual NH3 Clumps 



We ne xt look at the non- thermal line widths in the individual NH3 clumps in all three cores. 
In general, iJijina et al.l (|l999l ) found that non-thermal NH3 line widths in clustered environments 
are larger than those found in isolated regions. In Figure [TOl we plot Av^t versus Avt for the NII3 
clumps in each core. We omit NH3 clumps B3-A1 and C-Al wher e the correcte d non-thermal line 
width is effectively zero. We also show the best fit lines found by iJiiina et al.l to the relationship 
between thermal and non-thermal line widths in clustered and in isolated regions. Most of the Oph 
B clumps lie above the Avnt — Avt trend for isolated clumps and near the trend for the clumps 
in clustered re gions. The g ood agreement is somewhat surprising given that the majority of the 
objects in the Ijijina et al.l sample were observed with ~ 4 — 8 times poorer angular resolution, 
while those observed with high angular resolution are high mass star forming regions ~ 3 — 7.5 kpc 
distant and thus have very low linear resolution (excepting Orion B, at a distance of 420 pc). Even 
at the small spatial scales probed by our observations, Oph B is characterized by wide line widths 
that follow the relationship found for larger objects in clustered environments. In comparison, the 
NII3 clumps in Oph C lie well below the clustered Av^t — Avt trend, with Avnt < Avt, and also 
below the Avnt — Avt trend seen for objects not associated with a cluster. Two of the three Oph 
F clumps also fall significantly below the isolated object trend, while the third is more turbulent. 
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5.2. Discussion of the Cores 

5.2.1. Trends with N{H2) 

In Figure dH we plot the distribution of Tk, (Jnt / Cs, iV(NH3) and ^(NHs) with A^(H2) 
(calculated in §4.4 assuming = Tk) in Oph B, C and F. We additionally analyse Oph Bl and B2 
separately to examine any potential differences between the two. In Figure [TTk . we find that Tk 
values in Oph C are nearly universally lower than those in the other filaments, and show a tendency 
to decrease with increasing H2 column density. As described previously, the other Cores are warmer 
but also do not show a significant trend with A^(Il2). We show in Figure [TTb that Oph Bl and 
B2 are both consistent with having a constant, mildly supersonic ratio of non-thermal to thermal 
line widths over all A^(H2). Oph C line widths are generally dominated by thermal motions, and 
(^NT I Cs decreases significantly at A^(Il2) > 4 x lO^^cm"^. It is interesting to note that there are 
few data points at these column densities in the other Cores, and no observed decrease in a^T / Cs- 
Figure [TTb and d show that in all Cores A^(NH3) tends to increase with A^(Il2), following the 
general trend that the NII3 emission follows the continuum emission in the Cores. The differences 
between NH3 and continuum clumps are due to small scale differences in maxima. In addition, the 
fractional NII3 abundance, ^(NHs), tends to decrease with increasing A^(H2), although we note 
that Oph Bl and F have relatively few data points. Figure [TTb also illustrates the high A^(NIl3) 
values in Oph Bl relative to A^(Il2) compared with values in Oph B2, C and F. 

If NII3 is depleting at high densities, as the variation of X(NH3) with A^(Il2) suggests, NH3 
may not trace well the Tk in the densest and likely coldest regions. Hence, our assumption of 
Td = Tk may not be valid towards the highest columns and we may be underestimating and 



A^(H2) by factors of ~ 2. For example, IStamatellos et al.l (|2007l ) predicted that T^ in the centers 



of the Oph Cores could be as low as ~ 7K. Furthermore, X(NIl3) (oc A^(Il2), see Figure [TTI) may 
be overestimated by similar factors. High resolution multiwavelength continuum observations and 
models are needed to obtain independent assessments of T^ throughout the Oph Cores. 



5.2.2. aNT I Cs in Oph B 



Jijina et al.l (|l999l ) compiled numerous observations of NH3 in dense gas and found that while 
most starless clumps are characterized by largely thermal motions, a large fraction of clumps in 
clusters have (Jnt > and that the identification of an NH3 clump as part of a cluster has a 
larger impact on the observed line widths than association with a protostar. These findings are in 
agreement with our results in Oph B. The line widths in Oph B, with a mean Av = 0.89 km s~^ 
and (Ttvt / Cs = 1-5, are significantly wider and more dominated by non-thermal motions than those 
found in isolated cores. While B2 is associated with at least two embedded protostars and Bl 
appears starless, when studied separately (shown in Figure [T2]l . both Cores have similar aNT / Cs 
distributions. 
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If the non-thermal component is caused by turbulent motions in the gas, then it is interesting 
to consider the turbulence source. In the following, we consider the source of wide lines in this 
region as being due to "primordial" (i.e., undamped) turbulence, protostar-driven turbulence, bulk 
motions or biased sampling. 

Firstly, Oph B may have non-thermal motions throughout all the Core that are inherited from 
the surrounding cloud and that have not yet been damped. Since Oph B is associated with a few 
embedded protostars, however, it is likely that parts of the Core, at the very least, have been at high 
density for over a free-fall time tjf (~ 3 x 10^ years for n(H2) ~ 10^ cm~^, ~ nrr for the NH3 (1,1) 



and (2,2) transitions). Since the dissipation timescale for turbulence is ~ tff (jMac Low fc Klessen 



20041 ). it is unlikely that the non-thermal motions from the parent cloud have been retained, if 
the embedded protostars are indeed physically connected to the Core. Determining the relative 
velocities of the YSOs compared with the Core vlsr would help to address this question. 

Secondly, the embedded protostars in Oph B may be adding turbulence to the core through 
energy input associated with mass loss. One outflow has been found in CO emission in the re- 
gion associated with one of the two protostars near the peak NH3 integrated intensity in B2 
(IRS45/GY273 and IRS47/GY279, labelled in Figure H with a blue lobe towards the west and 



a red lobe towards the south-east (iKamazaki et al.l 120031 ). We do not, however, see localized re- 



gions of wide line widths associated with these or any protostars in Oph B. In fact, protostars are 
associated with Av minima. Additionally, wide line widths are found in Bl, where there are no 
embedded protostars and no known outflows. This suggests that the large non-thermal motions 
across the core are not driven by embedded YSOs. 

Thirdly, the wide line widths seen in the NII3 emission may be indicative of global infall 
in the Oph B Core. Given the mean NH3 line width in Oph B2, we calculate a virial mass 
Myir ~ 8M (Ti, assuming a dens ity distribution which varies as p oc r~'^. Given mass estimates of 



the Cores by lMotte et al.l (|l998l ) (which are uncertain to factors ~ 2) and accounting for the different 
Ophiuchus distance used by the authors (160 pc compared with our preferred value of 120 pc), we 
find M/Myir ~ 5 in Oph B2. Since the mean NH3 (1,1) opacity r = 1.6 over Oph B, the average 
individual hyperfine component is optically thin, and we would not expect t o see the asyrnmetr ically 



blue line profiles found in collapsing cores in optically thick line tracers. lAndre et al.l (|2007l ) find 
these spectroscopic signatures of infall motions in B2, with a clear blue infall profile towards our 
B2-A10 NH3 clump, and profiles suggestive of infall towards B2-A7 and other continuum objects 
in its eastern half. No evidence of infall motions in the tracers used were found towards central 
B2, but the lines they used (CS, II2CO and IICO+) may suffer from depletion at the high densities 
and low temperatures found at this location, masking any infall signature. Oph Bl, however, is 
also characterized by wide NH3 line widths and has M/Mmr ~ 1- Ad ditionally, spe ctroscopic infall 



signatures were found in Oph C towards our C-A3 NII3 clump by lAndre et al.l . where we find 
extremely narrow NII3 line widths (similarly, B2-A7 contains the narrowest lines observed in B2). 
It does not appear, then, that infall motions are the sole contribution to the wide lines we find in 
Oph B. 
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Lastly, rather than acting as a tracer of the densest gas in Oph B, the NH3 emission may be 
dominated in this high density environment by emission from the more turbulent outer envelope 
of the Core. Modelling of the expected emission given a constant abundance of NH3 would be 
necessary to determine accurately whether depletion is a factor in Oph B (note that we see evidence 
for decreasing abundance of NH3 with A^(H2); see Figure [TTll). Alternatively, observations of 
species which are excited at higher densities than lO^cm"'^, such as N2H"'" and deuterated species 
such as N2D"'" and H2D"'", may better probe the dense clump gas. Moderate resolution (~ 26") 
resolution observations of N2H"'" (1-0) (lAndre et al.ll2007l ) find narrower line emission in small scale 
N2H+ condensations (which we call 'clumps') in Oph Bl and the filament connecting Bl and B2 
{a NT = 0.15 lb 0.04 km s~^, or a^T / Cs 0.65 for Tk = 15 K), but non-thermal line widths in 
Oph B2 condensations remain transsonic (cjatt / ~ 1.1 < 2 on average). In an upcoming paper, 
we present H2D"'" lio — In observations in B2, showing that non-thermal line widths of gas at 
densities of ^ 10^ cm^^ (approximately the critical density of H2D'*", depending on the collisional 
cross section used) in the core are also transsonic, cjjvt ~ 0.25 kms~^, or a^T / Cs ^ 1.4 at 10 K (or 
~ 1.1 at 15 K; Priesen et al. 2009, in preparation). Regardless of whether or not the line widths 
seen in NH3 emission are tracing the highest density gas, the mean non-thermal motions of gas in 
B2 are greater than typically found in isolated cores. 



5.2.3. GNT / Cs in Oph C 



Motions in the Oph C Core differ substantially from those found in Oph B. Oph C is likely 
starless, like Oph Bl. In contrast to Bl (and also Oph B2 and F), Oph C is dominated by nearly 
thermal motions {{(Tnt I Cg) = 0.6), shown in Figure [T2j Oph C thus appears less affected by the 
clustered environment, and is kinematically more alike isolated clumps. NII3 clumps in the Pipe 
nebula and Perseus molecul ar cloud have similarly narrow line emission, for example, with the 
cnt / Cs ~ 1 — 2 in the Pipe (IRathborne et al.ll2008l 'l and intrinsic line widths ay typically less than 
0.2kms~^ in Perseus (jRosolowsky et al.ll2008al ). 



In addition, we see in Oph C a progressive decrease in the magnitude of non-thermal motions 
as we look at molecular lines that trace increasingly high densities. Non-thermal line widths 
o'NT = 0.20 — 0.28 k m s~^ were determ ined from DCO^ (1"0) observations (albeit with a large 



1'.5 beam) of Oph C (iLoren et al.lll99d ). At their peak emission locations in C, we find a^T = 
0.16 kms~^ and 0.17kms~^ for C2S and HC5N, respectively, from our GBT-only data. These 
locations coincide with the bulk of the NH3 emission in C, where we find the smallest a^T values 
(consistent at our velocity resolution wi t h nea rly or purely thermal motions), and where recent 
N2H+(l-0) observations by lAndrietay ([2OO3) also find aNT ~ 0.13 kms"^ (their C-MM3 - C- 
MM6). This trend indicates that turbulent motions have decreased at higher densities. 
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5.2.4- Oph B3 



Oph B3 is an unusual object within Oph B. It was originally detected in DCO^ emission 



(jLoren et al.lll99Cll ). It is not readily appar ent in dust co r itinuu m maps of Oph B, and consequently 
is only identified as a separate clump by IStanke et alj (|2006l ) . This lack of prominence suggests 
that interesting pot ential sites of star f o rmati on may be overlooked by purely continuum surveys. 
(In a similar result. iDi Francesco et alJ (|2004l ) found a region of extremely narrow N2H"'"(l-0) line 
width in Oph A associated with extended thermal continuum emission but not correlated with a 
continuum clump.) The relatively large difference in line of sight velocity between B3 and the Oph 
B mean (e.g., greater than that between Oph B and F) suggests that B3 may not be physically 
associated with the rest of Oph B. 

B3 is characterized in NH3 emission by moderately bright, extremely narrow lines which are 
nearly thermal in width. The lack of continuum emission implies a low H2 column density, leading 
to an NH3 fractional abundance lower limit similar to the X(NH3) values found for clumps in Bl 
and C. Interestingly, no C2S or HC5N emission was detected in B3, suggesting that B3 may be in 
a later evolutionary state than Bl and C. Alternatively, if B3 is physically distinct from the rest of 
Oph B, its initial chemistry may have differed. 

The kinetic temperature associated with the B3-A1 NH3 clump is 13.9 K. This is lower than 
the average Tk in Oph B but not particularly cold compared with other NH3 clumps in Oph B. 
The uncertainty in this value is large due to the low signal-to-noise ratio of the NH3 (2,2) line in 
the combined data, but we find a slightly lower value, Tk = 12.0 it 0.3 K, from the single-dish data 
alone. This temperature difference is not large enough to account for the lack of submillimeter 
emission at the B3 peak if the column density of material was similar to that found in Bl or B2. 



Based on the flux observed by IStanke et al.l (|2006l ) at the location of B3 (their MMS-108) 
and the gas temperature derived from our observations, we find a total clump mass for Oph B3 of 
M ~ 0.4 Mq, which is within a factor of ~ 2 (i.e., within uncertainties) of the virial mass calculated 
using the observed NH3 line width and core radius from Table [5] (M ~ 0.2 Mq). Accordingly, if 
it collapses, Oph B3 may form a very low mass star or brown dwarf, depending on how much 
material passes into a compact protostar. Future sensitive large format millimetre array detectors, 
such as SCUBA-2, should easily detect objects like B3 at greater than 3-cr levels of confidence with 
relatively short integration times. 



5.3. Implications for Clustered Star Formation 

Comparisons of the mean values of many parameters found in each of the Cores (Tables [6] 
and [7|) show significant differences between Oph B and C. In particular, Oph C is characterized 
by significantly narrower line widths and lower kinetic temperatures. This distinction is further 
illustrated in Figures [H] and US An object like Oph C (and perhaps Oph B3) resembles isolated 
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clumps in the close correspondence between the NH3 Une map and the dust continuum map, and 
in the evidence that the peak of the intensity map coincides with local minima in both Tk and At;. 
Such a core would fit well the idea that clustered star formation is just a spatially concentrated 
version of isolated star formation, with smaller, denser star forming clumps packed closer together 
than in isolated regions. The Oph B Core shows a remarkably different behavior. Oph B contains 
a larger number of NH3 clumps than Oph C or F, and there is much less correspondence between 
the locations of the NH3 and continuum clumps. The NH3 clumps show modest contrast with 
interclump gas in their intensity, and little contrast in their velocity or their line width. The 
greater amount of fragmentation of Oph B than found in C or F may be related to the relatively 
higher levels of turbulence in the Core, which are significantly greater than typically found in 
isolated regions. Thus, dense NH3 gas in Oph B does not resemble the dense gas in regions of 
isolated star formation, and this raises the issue whether this presents a different 'initial condition' 
for clustered star formation. 



6. Summary 

We have presented combined single-dish and interferometer NH3 (1,1) and (2,2) observations 
of the B, C and F Cores in the clustered star-forming Ophiuchus molecular cloud. We additionally 
present single-dish C2S (2i — Iq) and HC5N (9-8) observations of the Cores. Our main results can 
be summarized thus: 

1. While the large-scale features of submillimeter continuum emission and NH3 (1,1) integrated 
intensity appear similar, on 15" scales we find significant discrepancies between the dense gas tracers 
in Oph B, but good correspondence in Oph C and F. We find poor correspondence in Oph B between 
continuum clumps and NH3 clumps identified with 3D CLUMPFIND, with only five of fifteen NH3 
clumps located within 30" (3600 AU) of a dust clump. This is in contrast with previous findings of 
extremely high spatial correlation between the two gas tracers in isolated, low-mass starless clumps. 

2. We find vlsr varies little across any of the Oph Cores, and additionally varies by only 
~ 1.5kms~^ between them. 

3. Overall, the observed NH3 line widths in the Cores are generally large, and often slightly 
supersonic. We also find regions of localized narrow line emission (Av < 0.4 km s^^), some of which 
are associated with NH3 clumps. The larger line widths in Oph B {{Av) = 0.8 km s^^) agree with 
previous findings for clumps in clustered regions. Line widths in Oph C, however, decrease to nearly 
thermal values which are more representative of typical isolated clumps. 

4. The derived kinetic temperatures of the gas are remarkably constant across Oph B. Kinetic 
gas temperatures in both B and F are warmer {{Tk) = 15 K) than generally found in isolated star 
forming clumps, but are consistent with temperatures determined for cores in clustered environ- 
ments. The center of Oph C shows a minimum Tk ~ 9K, similar to previous results in isolated 
clumps. 
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5. We find no significant difference in Tk between NH3 clumps, submillimeter clumps ((TV) = 
14 K) and protostars {{Tk) = 15.5 K). Most other physical parameters have similarly insignificant 
variations between the objects. Wc do find that protostars are associated with significantly smaller 
line widths ((Au) = 0.4kms~^) with approximately equal contribution from thermal and non- 
thermal motions. 

6. We have determined NH3 abundance values towards the Cores, and find they agree with 
previous estimates of ^(NHs) in cold, dense environments. Single-dish observations of C2S and 
HC5N resulted in only a few detections, and derived column densities (~ 10^^ — 10^^ cm~^) similar 
to those found in other molecular clouds. 

7. It is unlikely that the wide line widths observed in Oph B are due to pervasive turbulent 
motions inherited from the parent cloud if the already formed embedded protostars are physically 
connected with the core. We find no evidence of influence by the protostars on the gas motions, 
i.e., through local increases in line widths or gas temperatures. The NH3 abundance in Oph B2 
appears to decrease with increasing gas density. We therefore suggest that the NH3 emission is 
biased, i.e., due to depletion at high inner densities, and is therefore not tracing the densest gas in 
the Core. This may explain the differences between the locations of NH3 and continuum clumps. 
Consequently, the gas temperatures may be lower in the center of Oph B than given by the NH3 
line ratios. 
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A. Determining Physical Parameters from HFS Line Fitting Results 

A.l. Kinetic Temperature 

Since the metastable states across i^-ladders are coupled only by collisions, and if we neglect 
the upper (J 7^ K) non-metastable states, then the populations of molecules in the metastable 
states can be described by the Boltzmann equation. 

n{J',K') _g{J',K') ( AE{J' , K'; J, K) \ 



n{J,K) g{J,K) "V Trot{J' , K'; J, K) 

where Trot{J' , K'; J, K) is the rotational temperature relating the populations in the (J, K) and 
(J', K') states, and AE{J', K'; J, K) is the energy difference between the two states. For (J', K') = 
(2,2) and {J,K) = (1,1), AE{2,2;l,l)/k = -41.5 K. The state statistical weights, q(J,K) an d 



g{J' , K'), are equal to 3 and 5 respectively for the NH3 (1,1) and (2,2) states (IHo fc Towneslll983l ). 
Assuming the molecular cloud is homogenous such that n{J',K') / n{J,K) = N{J',K') / N{J,K), 
we can solve for Tj-otC^, 2; 1, 1) using 



'^rot — TrotC^, 2; 1, 1) 



In 



-41.5[lnf- 0-^83 M2, 2) 
V r(l,l,m) Ai.(l,l) 

ATA{2,2,m) 
ATA{l,l,m) 



(A2) 



1 — exp(— t(1, 1, m)) 



n -1 



where AT^(1, 1, m) and Ar^(2, 2, m) are the antenna temperatures of the main component of the 
(1,1) and (2,2) transitions, respectively, T(l,l,m) is the opacity of the (1,1) line summed over the 
main component only, and we assume that the line widths Az^(l, 1) = Aiy{2, 2). We then calculate 
th e gas kinetic t e mper ature Tx from the rotational temperature following the updated result given 
by iTafalla et al.l (|2002l ) in their Appendix B: 



Tk 



rp21 

-'- rot 



%ln[l + l.lexp(-16/r2„i) 



(A3) 



A. 2. Column Density 



We calculate the column density of the (1,1) transition following iRosolowskv et al.l (|2008al ): 



iV(l,l) 



Svrt/g gl 1 
c2 52^(1,1) 



1 — exp 



kT 



(A4) 
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where gl = g2 are the statistical weights for the upper and lower states of th e (1,1) inversion 
transition. The Einstein A coefficient A{1, 1) = 1.68x 10"'^ s"^ (jPickett et al.lll998l ). and / T(z^)(ii/ = 
V2Tray/{cuQ)Ttot- The total NH3 column density A'^(NH3) can then be determined by calculating 
the value of the partition function Z of the metastable states: 



(A5) 



where the total NH3 column is then A^(l, 1) x Z/Z(l,l) . The values of the rotational constants 
B and C are 298117 MHz and 186726 MHz, respectively dPickett et alJll998l ^. The function S{J) 
accounts for the extra statistical weight of the ortho- over para-NHs states, with Z = 2 for J = 
3, 6, 9, ... and Z = 1 for all other J. 



B. Consequences of 0.3 km s ^ Velocity Resolution 

To fit the hyperfine components of the NH3 (1,1) inversion emission and obtain robust measure- 
ments of the opacity and excitation temperature of the line, we required an observational bandwidth 
wide enough to contain the main and at least two satellite components of the line. To meet this 
requirement, we were able to obtain velocity resolution across the band of only 0.3kms~^ due to 
the present correlator capabilities at the VLA. Here, we investigate the effect of the relatively low 
spectral line resolution on our determined line widths and kinetic temperatures by creating model 
spectra, including random noise with an rms value equal to that in our observed spectra. The 
spectra are then convolved by a finite spectral resolution and resampled. We then fit the resulting 
spectra with our custom HFS routine. 

NH3 (1,1) spectra were modelled given an excitation temperature T^x and opacity r typical of 
the regions observed, while the line width Av was varied. The line intensity at each velocity is then 
given by Equations [1] and [2j The amplitudes of the 21 components of the NH3 (2,2) line were then 



determ ined in a similar manner, with the relative intensities of each component given bv iKukolich 



(|l967l ). The relative intensity of the (2,2) line at a given Tk was determined using Equation IA2 



by first calculating the associated rotational temperature Trot (reversing Equation IA3j) . 

While the observational resolution does significantly impact some of the returned parameters, 
we show that the derived kinetic temperatures in this work are robust. 



B.l. Line widths 

Figure [13] plots the returned Av against that of the model spectrum for 0.1 km s^^ and 
0.3kms~^ resolution for a range of line widths. For both resolutions, the returned line width 
is greater than the actual for small values, with the 0.3kms~^ resolution showing a significantly 
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greater effect. The trends, however, fohow that expected for addition in quadrature of the true Hne 
width and the resolution, i.e. Avobs = V Avune + Avres with some additional offset at small values 
and larger scatter due to poor resolution. For the mean values found in this study, the corrected 
line widths consequently reflect the true line widths, but at small widths the corrected values still 
overestimate the true value. For this reason, we cannot accurately state the true widths of lines 
where returned (uncorrected) values Av < 0.35 km s~^, but do show that these uncorrected line 
widths are consistent with purely or nearly thermal values. 

B.2. Opacity 

Figure [HI plots the returned opacity r against that of the model spectrum for 0.1 kms~^ and 
0.3kms~^ resolution over a range of line widths. The relatively low spectral resolution of our 
observations has a significant impact on the returned opacity, with our fitting routine significantly 
underestimating the opacity by greater than 20% for line widths < 0.7 km s^^. In addition, the 
scatter in the returned opacities is large at small line widths, and this affects the uncertainty in the 
calculated NH3 column density, A^(NH3), in §4.3. Since A^(NH3) depends linearly on the opacity 
(see Equation IA4p , the uncertainty in the opacity dominates the uncertainty in the column density 
at small line widths. Since the column density also depends linearly on line width, however, and 
the returned line width at small values is greater than the true value by a similar relative amount, 
the returned column density is more accurate than the uncertainties suggest. 

B.3. Kinetic Temperature 

The results from the above tests show that at small NH3 line widths, the returned line widths 
and opacities can vary significantly from the true values. Regardless, the returned kinetic temper- 
atures from our NH3 line fitting are robust, as shown in Figure [15] (for Tk = 15 K). Even at small 
line widths {Av < 0.3kms~^), the kinetic temperatures derived from the fits are accurate within 
uncertainties of the true T^-, with relatively little scatter. At the lowest temperatures found in the 
Oph cores, Tk = 10 K, the results are similar with slightly more scatter at the lowest line widths. 
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Fig. 1. — The central region of th e Ophiuchus mol e cular c loud in 850 /im continuum emission 
originally mapped at the JCMT by Ijohnstone et alJ ( 2000bl ). Colour scale units are Jybeam"^, 
where the beam FWHM « 15". The B, C and F Cores are labelled (Oph B3 is not well detected 
in 850 /im continuum). Rectangles show areas mapped in NH3 (1,1) and (2,2) emission, as well 
as C2S (2i — lo) and HC5 N (9 - 8) at the GB T. Stars indicate locations of protostars identified in 
the infrared with Spitzer (jEnoch et al.l 120081 ). VLA and ATCA pointings were placed to provide 
Nyquist-sampled mosaicing of the indicated regions, with multiple beam overlap in the areas of 
bright continuum emission. 
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Table 1. Rest frequencies of all observed spectral lines 



Molecule Transition Rest Frequency 


Source 




GHz 




C2S (2i - lo) 


23.3440330 


Pickett et al. ri998) 


NH3 (1, 1) 


23.694495 


Ho & Townes ('1988) 


NH3 (2,2) 


23.722633 


Ho & Townes fl983) 


HC5N (9 - 8) 


23.963888 


Mvers et al. f 19791 


Table 2. GBT Observation Details by Region 


Core 


Area Mapped 


rms 




arcmin x arcmin 


K (Tmb) 


Oph Bl 


3 X 3 / 4 X 5.5 ^ 


0.05 / 0.07 ^ 


Oph B2 


5x4 


0.08 


Oph B3 


3x3 


0.07 


Oph C 


3x3/4x4^ 


0.05 / 0.07 ^ 


Oph F 


3x3 


0.06 



'^First values indicate the original OTF 
map size. Second values are the final map 
sizes, increased due to observed extent of 
emission. 

''First rms value is calculated over the cen- 
tral region where all observations overlap. 
Second value is calculated over the entire 
map. 
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Table 3. ATCA and VLA Observation Details by Region 







ATCA 




VLA 




Core 


Overlap Extent 




rms 


N ^ 


rms 




arcmin x arcmin 




mJy beam"-*^ 




mJy beam""-*^ 


Oph Bl 


2x3 


3 


20 


7b 


13 


Oph B2 


5x4 


10 


20 


16'= 


12 


Oph C 


4x4 


7 


30 


7 


10 


OphF 


3x2 


3 


30 


5 


15 



Note. — The spectral resolution of the observations is 0.1 kms ^ (ATCA) 
and0.3kms-i (VLA). 

^Number of individual pointings in mosaic observations. 

^Pointings for Bl were aligned and spaced to provide continuous coverage 
between Bl and B2. 

Additional pointings added to the B2 core to include B3 in the overlap 
region. 
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Fig. 2. — Observations of the Oph B Core (including Bl, B2 and B3). At lower left in each panel, ovals show the respective 
resolutions of the data, a) Integrated NH3 (1,1) main component intensity obtained with the GBT, ATCA and VLA telescopes 
with a synthesized beam FWHM 10.5" x 8.5". Emission was integrated over spectral channels with intensities > ±2 a rms. 
The colour scale is in Kkms^^ (Tr ). Contours begin at 3Kkms~^ and increase by 1.5Kkms^^. In both plots, stars de- 
note locations of Class 1 protostars ([Enoch et al.l 120081 ) . while triangles indicate t he positions of NH.s clu mps as identified by 
CLUMPFIND. (b) Continuum emission at 850 ^m in Oph B at 15" resolution as in I Jorgensen et al.l (|2008l ). The colour scale is 
in mJy beam~^. Contours beg in at 50mJybeam~^ and increase by 100 m Jy beam~^ . Squar es show locations of submillimeter 



to 



Enoch et al.l (j2008l ) 



cl umps (IJorgeiisen et al.ll2008l ). We also label the millimeter clump MM8 (jMotte et al.lll998l ) and protostars identified previous 
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Table 4. GBT C2S 2i - Iq and HC5N 9-8 Peak Parameters 



ID 


RA 


Dec. 


FWHM 


Vlsr 


Av 


Tmb 


N 


X 




J2000 


J2000 


AU 






K 


10^2 


10-10 


BI-C2S 


16 27 13.1 


-24 30 50 


5100 


3.67(2) 


0.33(5) 


0.36(3) 


4.3(7) 


3.1(7) 


C-C2S 


16 26 59.8 


-24 34 13 


9200 


3.95(1) 


0.33(2) 


1.12(4) 


10.6(8) 


1.5(1) 


C-HC5N 


16 26 58.8 


-24 34 30 


5800 


3.94(1) 


0.35(2) 


0.68(3) 


4.0(3) 


0.46(4) 
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Fig. 3. — Observations of the Oph C Core. At lower left in each panel, ovals show the respective 
resolutions of the data, a) Integrated NH3 (1,1) main component intensity obtained with the 
GBT, ATCA and VLA telescopes, summed over spectral channels with intensities > 2a rms. 
The colour scale is in units of Kkms~^ O^b)- Contours begin at 3Kkms~^ {Tb) and increase 
by 1.5Kkms~^. In both plots, triangles indicate the positions of NII3 clumps as identified by 



CLUM PFIND. b) Continuum emission at 850 ;um in Oph C at 15" resolution as in I Jorgensen et al 



20081). Contours and colour range as in Figure [21 Squares show locations of submillimeter clumps 



Jorgensen et al.l 120081 ). 



- 38 - 



39 00 



30 



40 00 



o 
o 

° 30 



41 00 



30 



-24 42 00 




16 27 30 27 24 21 

a (2000.0) 



27 24 
a (2000.0) 



IRS42/GY252" 

_j , 1 \ L 

21 18 



Fig. 4. — Observations of the Oph F Core. At lower left in each panel, ovals show the respective 
resolutions of the data, a) Integrated intensity NH3 (1,1) main component intensity obtained with 
the GBT, ATCA and VLA telescopes, summed over spectral channels with intensities >2a rms. 
The colour scale is in units of Kkms^^ (^b)- Contours begin at 3Kkm s~^ (Tr) and incr ease by 
IKkms^^. In both plots, stars denote locations of Class I protostars ( Enoch et al.l 120081 ) . while 
triangles indicate the positions of NII3 clum ps as identified b y CLUM PFIND. b) Continuum emission 
at 850 /im in Oph F at 15" resolution as in lJorgensen et al.l ( 20081'). Contours and colour range as 
in Figure [2j Squares show locati ons of submil l imete r clumps (jJorgensen et al.ll2008l ). We also label 
protostars identified previous to Enoch et al.l (j2008l ). 
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Fig. 5. — a) Spectra of all species observed at the GBT at the integrated intensity C2S (2i — Iq) 
peak in Oph Bl. The NH3 (1,1) and (2,2) and C2S baselines are offset from for clarity, b) Spectra 
of all species observed at the GBT at the integrated intensity C2S (2i — Iq) peak in Oph C. The 
NH3 (1,1) and (2,2) and C2S baselines are offset for clarity. 



-40- 



0.1 0.2 0.3 0.4 0.5 




a (2000.0) a (2000.0) 



Fig. 6. — GBT observations of Oph C. The circle in each panel gives the relative FWHM of the 
GBT data, a) The Ophiuchus C core in 850 ;uni continuum emission at 15" resolution (greyscale, 
Jybeam"^ units) overlaid with continuum contours (grey). Black contours show the integrated 
NH3 (1,1) intensity at ~ 32" resolution, beginning at 3Kkms~-^ and increasing by 1.5Kkms~^. 
For all plots, emission was integrated over spectral channels with intensities > ±2(7 rms. b) 
850 //m continuum emission in Oph C overlaid with black NH3 (2,2) integrated emission contours. 
Contours begin at 0.3Kkms^^ and increase by O.lKkms"^. c) 850 /um continuum emission in 
Oph C overlaid with black C2S 2i — Iq integrated emission contours. Contours are 0.2, 0.25, 
0.3Kkms~^. d) 850 //m continuum emission in Oph C overlaid with black HC5N 9 — 8 integrated 
emission. Contours are 0.1, 0.2Kkms~^. 
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Table 5. NH3 (1,1) clumpfind peaks and parameters 



ID 


RA 


Dec. 


FWHM^ 


X 


FWHM^;^ 


Peak 






J2000 


(AU) 




Bl-Al 


16 27 9.7 


-24 


27 43.0 


3900 


X 


6700 


4.90 


B1-A2 


16 27 12.4 


-24 


29 49.0 


6800 


X 


7800 


5.10 


B1-A3 


16 27 15.4 


-24 


30 29.1 


5300 


X 


5000 


5.60 


B1-A4 


16 27 16.7 


-24 


30 59.1 


3700 


X 


2900 


4.90 


B2-A1 


16 27 13.8 


-24 


27 35.0 


2300 


X 


2500 


5.20 


B2-A2 


16 27 16.0 


-24 


27 39.0 


5600 


X 


6100 


5.10 


B2-A3 


16 27 19.0 


-24 


27 11.0 


4700 


X 


6700 


3.90 


B2-A4 


16 27 24.2 


-24 


27 29.0 


6400 


X 


6000 


5.10 


B2-A5 


16 27 25.0 


-24 


26 55.0 


5300 


X 


4700 


5.00 


B2-A6 


16 27 29.0 


-24 


26 59.0 


3800 


X 


6700 


5.70 


B2-A7 


16 27 32.5 


-24 


26 57.0 


3800 


X 


5200 


6.70 


B2-A8 


16 27 32.8 


-24 


26 17.0 


3600 


X 


3600 


4.00 


B2-A9 


16 27 33.5 


-24 


27 43.0 


5000 


X 


3600 


3.60 


B2-A10 


16 27 33.8 


-24 


26 14.9 


2700 


X 


3700 


4.00 


B3-A1 


16 27 21.8 


-24 


24 42.9 


3900 


X 


6100 


3.60 


C-Al 


16 26 59.0 


-24 


34 14.8 


5000 


X 


5600 


5.40 


C-A2 


16 26 59.1 


-24 


32 54.8 


4500 


X 


4400 


4.30 


C-A3 


16 27 1.2 


-24 


34 34.8 


6000 


X 


5400 


5.50 


F-Al 


16 27 21.9 


-24 


39 52.0 


3600 


X 


3200 


4.50 


F-A2 


16 27 24.2 


-24 


40 52.1 


4400 


X 


5900 


4.90 


F-A3 


16 27 27.1 


-24 


40 58.1 


2400 


X 


5000 


5.50 



'^FWHM values calculated by clumpfind have not been decon- 
volved with the synthesized beam of the combined data. 
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Table 6. NH3 (1,1) Line Characteristics in Combined Data 



Core 


Value 




Mean 


RMS 


Min 


Max 


OphB 


Visr (kms" 


-1-1 


3.96 


0.24 


3.15 


4.57 




Av (kms~ 


-1-1 


0.83 


0.21 


0.08 


1.37 




r 




1.6 


0.8 


0.5 


4.7 




Te. (K) 




9.9 


2.3 


5.4 


19.5 


Oph C 


Visr (kms" 


-1\ 


4.01 


0.07 


3.81 


4.10 




Av (kms~ 


-1\ 


0.37 


0.13 


0.11 


0.70 




r 




4.7 


2.6 


0.7 


10.8 




Te. (K) 




6.9 


1.2 


5.1 


11.6 


OphF 


Visr (kms" 


-1\ 


4.34 


0.12 


4.24 


4.82 




Av (kms~ 


-1\ 


0.53 


0.30 


0.10 


1.13 




T 




1.9 


0.8 


0.7 


3.6 




Tex (K) 




9.5 


2.5 


5.0 


14.3 
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Fig. 7. — a) Line velocity or vlsr in Oph B. Colour scale is in kms~^. The values shown are those 
derived after convolving the data to 15" resolution and then regridding into 15" x 15" pixels. 
In all plots, contours show integrated NH3 (1,1) intensity convolved to 15" resolution, beginning 
at 4.5Kkms~^ and increasing by 1.5Kkms~^. Stars indicate protostar positions and triangles 
indicate positions of NH3 clumps, b) Fitted Av in Oph B. Colour scale is in kms~^. c) Tk in Oph 
B. Greyscale from 10 K to 20 K. d) Ratio of the non-thermal to thermal line width components. 
Greyscale from aNT / = to 2.4. e) Total column density of H2 derived from 850 //m dust 
continuum observations in Oph B, regridded to match the combined NH3 observations. The dust 
temperature per pixel was assumed to be equal to the gas temperature derived from HFS line 
fitting of the NH3 observations. The A^(H2) values shown have been divided by 10^^. f) Fractional 
abundance of NH3 divided by 10~^. 
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Fig. 8. — a) Line velocity or vlsr in Oph C. Colour scale is in kms~^. The values shown are 
those derived after convolving the data to 15" resolution and then regridding into 15" x 15" 
pixels. In all plots, contours show integrated NH3 (1,1) intensity convolved to 15" resolution, 
beginning at 4.5Kkms~^ and increasing by 1.5Kkms^^. Stars indicate protostar positions and 
triangles indicate positions of NH3 clumps, b) Fitted Au in Oph C. Colour scale is in kms~^. 
c) Tfc in Oph C. Greyscale from 10 K to 20 K. d) Ratio of the non-thermal to thermal line width 
components. Greyscale from gj^t / Cs = to 2.4. e) Total column density of H2 derived from 850 //m 
dust continuum observations in Oph C, convolved to 15" resolution and regridded to match the 
combined NH3 observations. The dust temperature per pixel was assumed to be equal to the 
gas temperature Tk derived from HFS line fitting of the NH3 observations. The A^(H2) values 
shown have been divided by 10^^. f) Fractional abundance of NH3 divided by 10"^. 
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Fig. 9. — a) Line velocity or v^sR in Oph F. Colour scale is in kms^^. The values shown are those 
derived after convolving the data to 15" resolution and then regridding into 15" x 15" pixels. 
In all plots, contours show integrated NH3 (1,1) intensity convolved to 15" resolution, beginning 
at 3Kkms~^ and increasing by IKkms"^. Stars indicate protostar positions and triangles indi- 
cate positions of NH3 clumps, b) Fitted Av in Oph F. Colour scale is in kms~^. c) T^, in Oph 
C. Greyscale from 10 K to 20 K. d) Ratio of the non-thermal to thermal line width components. 
Greyscale from a^T / Cg = to 2.4. e) Total column density of H2 derived from 850 //m dust 
continuum observations in Oph F, convolved to 15" resolution and regridded to match the com- 
bined NH3 observations. The dust temperature Ta per pixel was assumed to be equal to the gas 
temperature derived from HFS line fitting of the NH3 observations. The A^(H2) values shown 
have been divided by 10^^. f) Fractional abundance of NH3 divided by 10^^. 
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Table 7. Physical Properties of Filaments Derived Prom Fitted Parameters 



Core 


Value 


Mean 


RMS 


Min 


Max 


yj^jiL j_> 




1 ^ 1 

i-KJ. 1. 




19 S 






aNT (kms~^) 


0.35 


0.09 


0.04 


0.57 




O'NT 1 Cs 


1.5 


0.4 


0.5 


2.5 




Ar(NH3) (x IQi^'^cm-'^) 


22 


8.8 


1.3 


48 




n(H2) (x 


5.9 


10 


0.8 


77 




X(NH3) (x 10-9) 


14 


9.1 


2.4 


41 


Oph C 


Tk (K) 


12.0 


1.6 


9.4 


19.8 




gnt (kms"^) 


0.14 


0.06 


0.0 


0.29 




CATT / Cs 


0.6 


0.2 


0.0 


1.2 




Ar(NH3) (x 10^3cj^-2) 


24 


11 


7.1 


55 




n(H2) (x lO-^cm-^) 


2.0 


1.0 


0.8 


4.2 




^(NHs) (x 10-9) 


8.2 


4.4 


2.8 


28 


OphF 


Tk (K) 


15.5 


2.5 


12.8 


22.9 




CWT (kms"^) 


0.21 


0.13 


0.05 


0.47 




catt / Cs 


0.9 


0.6 


0.2 


2.2 




Ar(NH3) (x 10^3cjn-2) 


14 


5.9 


3.5 


25 




n(H2) (x lO-^cm-^) 


3.6 


2.4 


0.9 


10.6 




X(NH3) (x 10-9) 


5.4 


1.8 


2.2 


10 



Table 8. 



Derived Parameters at NH3 (1,1) Peak Locations 



ID 


Vlsr 


Av 
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Tk 
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^ ex 


a NT 


7V(H2) 


iV(NH3)tot 


X(NH3) 


n(H2) 




kms~^ 






K 


K 


kms~^ 


10^^ cm-2 


10^^ cm-2 


10-9 


10"^ cm-'^ 


Bl-Al 


4.Ul(dj 


U.OO{l) 


l.D(3j 


lD.3(1.9j 


c\ /I / 1 ^ 
9.4(lj 


0.35(3j 




00 / /I ^ 
z2(4j 




2.8(1) 


"D 1 AO 

rSi-Az 


o.o4(^Z ) 


u.yz^o j 


1 ^ \ 
i.O(i j 


iD.4(i.O j 


12.1(1) 


U.oo(Z j 


14(3) 


on/ A \ 
Zy(4j 


17(8) 




JdI-Ao 


o.yo^i j 


n ^^o/'Q^ 
U.DZl^O j 


Q /I /'0\ 

o.4(Z j 


io. /(i.U j 


9.7(1) 


n ^ /oA 
U.Z5(Z j 


^ Q{ A\ 
18(4) 


41(o j 


OQ /'n^ 


/I 0/1 \ 

4.2(1) 




Q CQ ^0\ 

o.oo[Z ) 


U.D4(^o J 


/I 0/'Q\ 

4.Z(o ) 


10 /I /n ^^ 

iz.4(u.y J 


8.4(1) 


n oo/'o^ 
U.ZZ(Z j 


i2(2j 


4Z(0 j 




3.2(1) 


B2-A1 


4.18(2) 


0.57(5) 


1.6(3) 


13.3(1.2) 


11.6(1) 


0.23(2) 




22(4) 




13(1) 


"DO A 


A 10/ 0\ 

4.io(Zj 


0.8d(o j 


O.o( j 


14.1(=^) 


in /1 ^ 
iy.o(l j 


O.o5( j 




1»( ) 




0-0( ) 


T30 A 

Jdz-Ao 


/I 1 n/'Q\ 
4.iU(oj 


U.oU( r j 


l.a(aj 


"1/1 o^ 
i4.z(i.zj 


in /'i ^ 
l(J.z(i j 


O.oo(oj 


15(2) 


on/ /I ^ 
z0(4j 


14(6) 


4.7(1) 


T30 A A 


A nn/'o\ 
4.UU(zj 


l.UU(o j 


1 /I /oA 
i.4(z j 


lo. 5(1.0 j 


11 /' 1 ^ 
11.6(1) 


n /1 /'Q^ 
0.4z(dj 




QO/ 

oz(5j 


C\ ^ I A\ 

9. r (4j 


1 1 'y { q\ 

11.7(8) 


"DO A K 

JdZ-Ao 


Q ni /' oA 
o.yi^Z j 


u. ^y(o j 


/I /oA 

z.4(z j 


1 Q 0/ 1 n\ 
io.Z(i.U j 


9.7(1) 


n Qo /oA 
U.oZ(Z j 




QQ / /I \ 

oo(4j 


f;/Q\ 
O.D(0 J 


4.8(2) 


"DO A 

Jdz-ad 


4.oZ(^Z j 


U.Db(^4J 


1 Q /'o\ 
i.o(Z j 


14.5(1. i ) 


1/1 1 /I ^ 
14.1(1) 


n oT/ o\ 
U.z ( (zj 


33(7) 


0/( //( \ 

z4(4 j 


D.D(oj 




"DO A T 


4.57(1 ) 


n oo/'o^ 
O.oo(z j 


00/' o\ 
z.o(zj 


13.9(1.1 j 


1 1 n/1 ^ 
11.9(1} 


n 1 o/'o\ 
0.1z(zj 


35(7) 


oo/o\ 
zz(zj 


fi K /o^ 
b.5(z j 


1 1 tifQ\ 


"DO A 


A OKfO\ 

4.oo(oj 


0.95(8; 


0.5( j 


14.D( j 


16.9(1 j 


O.o9( j 


43(7) 


If?/ *\ 

lo( j 


3.7(=^) 


0-0( ) 


B2-A9 






















B2-A10 


4.22(2) 


0.71(6) 


1.6(3) 


14.3(1.2) 


10.3(1) 


0.29(3) 


35(7) 


23(4) 


6.1(3) 


4.8(2) 


B3-A1 


3.15(2) 


0.08(8) 


0.7(*) 


13.9(*) 


13.9(1) 


0.0(*) 




1.3(*) 




12.4(*) 


C-Al 


4.08(1) 


0.16(1) 


9.3(6) 


10.4(0.7) 


8.0(2) 


0.0(*) 


84(14) 


33(2) 


4.0(2) 


4.2(2) 


C-A2 


3.97(1) 


0.36(4) 


2.7(6) 


12.0(1.2) 


8.8(2) 


0.13(1) 


14(3) 


20(4) 


14(9) 


4.1(2) 


C-A3 


3.98(2) 


0.29(2) 


9.3(5) 


10.7(0.8) 


8.0(2) 


0.10(3) 


68(4) 


60(5) 


8.8(2) 


3.9(2) 


F-Al 


4.82(4) 


0.8(1) 


1.5(7) 


15.3(2.0) 


7.4(2) 


0.32(5) 


31(6) 


16(7) 


5.1(4) 


1.6(1) 


F-A2 


4.31(2) 


0.58(5) 


1.8(4) 


15.6(1.6) 


9.9(2) 


0.23(2) 


35(7) 


19(5) 


5.5(3) 


3.4(1) 


F-A3 


4.31(2) 


0.43(3) 


2.6(4) 


16.8(1.8) 


9.8(2) 


0.16(2) 


24(5) 


19(3) 


7.8(3) 


3.0(1) 



Note. — Uncertainties are given in brackets beside values and show the uncertainty in the last digit of the value, 
with the exception of the Tk uncertainties which are given in K units. A (*) indicates that the uncertainty in the 
value at this position were large, even though in most cases the values themselves are reasonable. A good HFS fit 
was not found at the peak position of NH3 clump B2-A9. 
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Fig. 10. — Non-thermal versus thermal line widths (FWHM) for individu al NH3 clumps in the 
Oph B, C and F Cores. The two dashed lines show relationships found bv iJiiina et al.l (|l999l ) for 
NH3 clumps both with and without an associated cluster. Clumps in Oph B follow the clustered 
trend, but clumps in Oph C are nearer the relationship for more isolated clumps. Clumps in Oph 
F scatter across the plot. The dotted line represents Avnt = Avt- 
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Fig. 11. — a) Kinetic temperature (Tk) in K versus A^(H2) (derived from the 850 ^m submillimeter 
continuum data) in cm~^ units for Oph Bl, B2, C and F. In all plots, points shown are values in 
15" pixels. Temperatures in Oph C are nearly universally lower than those in Bl, B2 and F, and 
show a tendency to decrease with increasing H2 column density. A significant trend in temperature 
with A^(H2) is not obvious in the other Cores, b) a^T / Cg versus A^(H2) in Oph Bl, B2, C and 
F. The dashed line indicates a^T /cs = 1. Oph Bl and B2 are consistent with a constant ratio 
of non-thermal to thermal line widths, where a^T / Cs > 1, over all H2 column densities. Oph C 
line widths are generally dominated by thermal motions, and ajsfT / Cg decreases significantly at H2 
column densities above log(A^(H2)) ~ 22.6. c) A^(NH3) (derived from the NH3 HFS line fitting 
results) versus A^(H2) for Oph Bl, B2, C and F. NH3 column densities tend to increase with H2 
column densities in all Cores. Oph Bl has significantly higher A^(NH3) for its relatively low A^(H2) 
values compared with the other Cores, d) X(NH3) versus A'^(H2) for Oph Bl, B2, C and F. While 
c) shows the NH3 column densities tend to increase with increasing A^(H2), the fractional NH3 
abundances appear to decrease with higher H2 column densities. 
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Fig. 12. — Histogram crjvr / Cg in Oph Bl, B2, C and F. Values plotted are those in 15" pixels. 
Oph C is characterized by significantly lower non-thermal to thermal line width ratios than the 
other Cores. When analysed separately, Oph Bl and B2 are very similar in a^T / Cs- 
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Fig. 13. — Input versus recovered line widths determined by creating a model spectrum convolved to 
0.1 km s^^ and 0.3 km s~^ velocity resolution and subsequently fitting with the NH3 HFS routine. 
Solid lines show the expected observed line width Avobs = V ^vune + Avres for O.lkms"^ and 
0.3 km s-^ 
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Fig. 14. — Recovered versus input NH3 (1,1) opacities as a function of line width Av for 0.1 kms"-*^ 
and 0.3km s~^ velocity resolution. At small line widths, the returned opacities are significantly 
smaller than the true values. 
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Fig. 15. — Recovered versus input kinetic temperatures Tk as a function of NH3 (1,1) line width 
for O.lkms"^ and 0.3 km s""^ velocity resolution. Even at small line widths, the returned kinetic 
temperature closely agrees with the input value, with little scatter. 
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